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Este número de la Revista SELPER aparece con un 
retraso significativo, producto de las. dificiles cir 
lancias económicas y polílicas de nuestra 
américa, slluación que necosariamónte ro- 
percuta en nuestra Sociedad. 
Sin embargo, debido a un gran esfuerzo de cont 
nuidad realizado por las Sede SELPER Perú y la Sede SELPER Colombia, coordinado con el 
Comité Editorial dde SELPER, ha sido posible finalmente producir esta edición, quo creamos será 
do gran interás y ullidad para nuestra comunidad científica, dada la calidad de los artículos 





NOTA DEL EDITOR 
Prof: MAURICIO ARAYA-FIGUEROA 
Presidente SELPER 1983 1986. 














incluidos, 
Destacan las contribuciones del Dr. Jorge Lira (UNAM, México) quien propone un Sistema de 
3 Digital de facil uso y que está disponible, on forma gratuita, para todos los Socios de 

uye un interesante Análisis Digital cdo Texturas en 

obrando una gran importancia, debido al lanzamion- 
Debido a la trascendencia de estas dos contribuciones, se 
ón en inglés, para dar una mayor amplitud do aplicación 





soliciten, También se li 
Imagenes de Radar (SAR); tema que est 
to del Satélite Europeo ERS-1. 
¡solicitó al Autor la inclusión de una ve 
a estos trabajos de investigación 
Sé presentan, además, dos notables contribuciones de Investigadores del CLIRSEN, Ecuador, 
le Imagenes LANDSAT/TM y ERS-1/SAR. Finalmen- 
>, que Se abre con esta nueva herramienta 

















sobre Erupciones Volcánicas y! 


1o, para completar e 


























de Imágonos de Radar, se incluye una breve reseña sobre resultados goneralos obtenidos hasta 
ahora con el sistema ERS-1. que representará un excelente complemento a la Información 
LANDSAT, SPO sistemas normalmente utilizados en Percepción Remota. 

Creemos que e antoner la continuidad de las Publicaciones SELPER, aunque 
existan retrasos por causas q a voluntad, El esfuerzo que realizan 


los Autores para contribuir con sus investigaciones, ia labor dol Comité Editorial Internacional 
por sallr adelante con estas adiciones, el importante apoyo de las Instítuciones Latinoamerica: 
nas (IGM, ONERN, IGAC, INPE y otras) e Intemacionales (ESA, DLA, CNES, MDA, F'S, 
ASHTECH, INTERGRAPH, INTERA, ITC y otras). bien merecen verso materíalizados an es! 

ediciones, que son el mejor ejemplo de -coperación internacional efectiva realizada 














durante tantos años 
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SISTEMA DE ANALISIS DIGITAL 
DE IMAGENES, V. 2.0 . 


Jorge Lira y Rosario Landgrave 


Insituto de Geofisica-UNAM - A. P. 22-560. 
14000 México D.F. 








Resumen 


La versión 1.0 del Sistema de Andiisis Digital de Imágenes (SAND!) anteriormente 
reportada (Lia tal, 1983) ha sido rediseñada y reeditada, La versión: 2.0 actual porta más 
de 70 módulos operados por menús. Todos los módulos están escritos en TurboPascal 6.0 
operando bajo DOS 5.0. Se requiere de 1 MB RAN, coprocesador matemático y ratón, 
además de un disco duro para almacenar las imágenes digitales. El sistema, trabajando en 
una PC 386 6 486 ostá dividido en dos grandes grupos: Procesamientos y Realces y 
Reconocimiento de Patrones. Un grupo de utlerías para desplegar y manipularlas imágenes. 
se encuentra disponible en varios lugares del paquete, La ventana de despliegue es 
desplazable empleando llaves preprogramadas conjuntamente con varias paletas de color 
incluyendo blanco y negro. Las imágenes son desplegadas a través de tarjetas de 
despliegue comerciales y son almacenadas en formato libre DOS, estando lógicamente 
relacionadas con etiquetas de control que incluyen información sobre el tamaño de la 
imagen, historia del procesamiento y un comentario del usuario. Las principales 
Aplicaciones están dedicadas ala Percepción Remota. Este sistema está disponible, bre de 
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costo, a la Comunidad de Especialistas, en particular ligados a SELPER. 





INTRODUCCION 


a tasa de alía productividad en nuevos 

métodos de procesamiento digital de 

Ímágenos y reconocimiento de patrones ha 

ido sustentada en a Úlima década: el mo 
GUA NO muestra signos do desinar, por el contara 
úfuévas revistas y aniculos especializados emergen de 
manera constanto. Do aquí que un mec ampio para 
probar nuevos modelos y mátados de anália de 
imágenes so vea deseable. Este medio debe 1raducirseen 
Un sistema do propósito general, de cardeter modular y la 
más transportablo posible. El gran desarro en e campo 
¿e las microcomputadoras, combmado con las nuevas 
tratos de despliegue comerciales de at resolución. 
sugiera que una PC-380 0 486 sea una platatorma 
adecuada para sopoñar tal ustema. El sistema dote ser 
Jaciment actuaizable ya que su penal objeivo es. 
apoyar las. actvidades de Investigación en el andiisis 
“gal de imágenes, Cada adición o actualización de un 
médula incrementa la capacidad de análisis del sistema, 
oyendo beneícios adicionales al proceso tales coma, 
proporcionar capacidad de aplicaciones rutinarias. Al 
imponer un. sistema basado en una PO se restringe 
aLomálicamente el uso de cies sistemas operatwos y 
defie un número de lenguajes a ser empleados. 

No obstanto, el progreso en velocidad, desanaño de 
iemori, compladores y sistemas operativos para PCs 
predice la disponblidad de un amplo mecio para el 
procesamiento de imágenes en un luuro proximo. 














Esta no es una sta completa como tampoco una recomendación 






DRAGON)", cuyos precios van desde mi 
decenas de mies de dolares y trabajando en varie 
plataformas y medios computacionales. Algunos oros 
paquetes están publicados en la iteratura con diferantos. 
capacidades y medios de operación (Borstein y 
Formeyhough, 1975; Wharton y Turner. 1981; Landy et 41 
1984; Harrington er al, 1986; Bryant, 1990, Sador y Winne, 
1991) incluyendo algoritmos para gráficas y 
procesamiento de imágenes (Pavkcis, 1982), Algunos do 
llos están ospecisizados hacia una tarea especifica 
(Warton y Tumer, 1981; Bryant. 1990) y otros hacia un 
¡campo deferido como la Percepción Remota (Bemstein y 
Femeyhough, 1975; Landy et al. 1984), un sistema 
expe (Schowengert y Wang, 1989), o para educación 
(Sader y Wine, 1991) 

Presentamos aquí un sistema modular de propósito 
bajo MS DOS 
5.0, fácilmente expandible y totalmento escrito en 
TurboPascal 6.0, cuyo Objetivo principal es apoyar las 
actividades de investigación y desarrollo en ol área de 
análisis digital de imágenes, aún cuando Lambién es 
empleado para actividades educacionales. Nuestro 
sistema presenta las siguientes caracteristicas: 1. Su 
diseñoes funcional, 2- Está estructuradoen cuantos tareas 
de realces y reconocimiento de palrones, 3 - Su operación 
es amigable y 4.- Es delácilactualización, 5-SANDIes una 
plataforma para soportar nuevos desarrollos, probar 
métodos y algoramos enel área de análisis de Imagenes en 
general. La estructura y operación de esto sistema está 
descrita en tas siguientes ds secciones. 
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Estructura de SAND! 


Antes de cuantificar las propiedades de un patrón es 
generalmente requerido el realzario o ponerlo en 
evidencia. En muchas ocasiones un patrón o una ciase de 
patrones no está bien diferenciado con respecto al resto 
de la imagen, ruido, borrosidad, interacción interpixel o 
inhomogoneidades en a luminación pueden oscurecer u 
ocultar la presencia de un patrón haciendo dificil su 
identificación, Un realce es en general requerido para 
poner an evidencia al patrón previo al proceso de 
reconocimiento; esta suposición básica es la que define la 
estructura (Fig. 1) de SANOL Los Procesamientos y 
Realces son clasificados de acuerdo al tipo de operación 
que emplean, así entonces los Roaices Radiométricos 
están para la manipulación del histograma de la imagen. 
Los Rualces Geométicos son para poner en evidencia 
algunos aspectos geométricos de la imagen: Los Realces 
de Fidelidadson básicamente restauraciones ala imagen, 
ya que tienden a aproximar, lo más ficimente posidle, una 
imagen degradada a un modelo ideal de la misma. 
Reconstrucción y Compresión utliza zoom por solnes e 
interpolación lineal y reducción de la imagen por 
'decimación, es decir alÍminando un pixel si y otro no: Las 
Operaciones al Conjunto son empleadas para transformar 
Una imagen multespeciral como un todo. Cada una de 
stas secciones Incluye un conjunto de operaciones de 
acuerdo al diagrama (Fig. 1) de SAND! 

La Sección de Reconocimiento de Patrones está 
subdividida en Patrones Espaciales y Espectrales, aún 
cuando entro estas secciones existe una relación 
estrecha, La segmentación de patrones espectrales 
incluyo una técnica por densidad del histograma. un 
clasificador no-superisado utilzando el Algorimo de 
Isodata(Touy Gonzalez, 1974) y un clasificador Bayesiano 
(Fichards, 1986) empleando el principio de máxima 
verosimilitud, Un Algoritmo de Crecimiento de Regiones 
(Wani, 1987) está incluido conjuntamento con los 
úsigulentos prodicados de unttormidad: a) - Diferencia de 
Promedio, b).- Correlación yc). Entropía, el primero para 
"segmentar portonos y los dos úlimos para segmentar por 
texturas. Un analizador contoxtual basacio én la matriz de 
co-9curencia y con el criterio de entropía trabaja para 
Imágenes monoespectrales y es especificamente 
adecuado para análisis de imágenes de radar. En lo que 
concierna a Patrones Espaciales, el paquete incluye 
segmentación orientada hacia la detección de lineas, 
manipulaciones algebraicas de bandas, operaciones para 
Imágenes clasiicacias donde cada clase se maneja como 
un conjunto y morfología matemática (Serra, 1982) en dos 
y los dimensiones, conjuntamente con aplicaciones 
selectas. Un grupo de Utlerias se incluye en lugares 
selectos del paquete para apoyar la manipulación de las. 
imágenes 

Internamente, SANO! está distribuido en programas y 
Unidades del TurboPascal, los menús forman un grupo 
aparte de programas. Las unidades agrupan un cierto 
número de procedimientos para aplicaciones 
especializadas. otras unidades incluyen procedimientos 
para manejo de archivos, reportes y manejo del 
despliegue de las imágenes, La adición de nuevos 
algormos involucra únicamente la actualización del menú. 
y la unidad donde el programa o procedimiento es 
Incluido la recompilación por anto está restringida a ésta. 
part del sistema. El control de menús, los reportes sobre 
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estores y el progreso de una operación es llevada cabo 
por procedimientos especialmente diseñados que 
pueden ser fácimente modificados para. acomodar 
nuevas tareas. E 

Las imágenes son almacenadas on formato ¡bre DOS 
(raster) con un tamaño de bloque Igual al número de 
columnas de la imagen: las etíquetas son manejadas en 
renos separados con bloques de 128 bytes: entre las 
etiquetas y las imágenes oxiste una relación lógica, El 
"máx tamaño de ls imágenes está resiigldo a 1024 
1024 pirels básicamente debido a las limitaciones que 
imponeolDOS Parala Transformada do Fourier, eltamaño 
está mado a imágenes de orden 2 (con n entero) y el 
crecimiento de regiones a. 400 x 400 pixelo. Los 
úiasilicadores especiales, supervisado y no-super- 
visado, operan con 7 bandas y 10 clases máximo, El 
analisis canónico, por otrolado,rabajahasta para Sciases 
y 7 bands. Para imágenes de 1024 x 1024 50 emplea un 

de Transpuesta Rápida (Huang, 1981) on la 
Transformada de Fourler, ya que ésta sa encuenta 
implantada como una cascada de dos transformaciones 
unidmensicnales. 

Se emplea un byto por pitel en las imágenes, con 
excepción de aqueñas temporales. donde, con varios 
grados de precisión pueden requerirse dos, cuatro ocho 
bytes por pixel, Los archivos temporales son 
automáticamente borrados una vez que la operación ha 
torminado. Las etiquetas de contol (Fig, 2) necesanas 
para identificar una imagen son generadas empleando 
una delas eras, Las atíquetas incluyen el tamaño dela 
imagen, número de pirels por Inga y número de Incas, 
número de procesamientos recibidas y un comentario por 
+ usuario de 80 caractores máximo. Un comentario de la 
histosa de procesamientos aplicados a la imagen os 
“agregado cada vez que una operación os aplicada ea 

En a siguiente gura es un entro y im número, 
¿e alquetas de historia y im (mod 2) = O, ortoncos Km 
(6) + 1, de ora manera K = md) +2. 





Lasimágenes son identificadas por sunombre y pora 
extensión POL (Procesamiento Digital de Imágenes). 
ú“aquelas sin etiquetas llevan la extension IMA (Imagen). 
Los archivos de etiquetas tienen el mismo nombre y 
“extensión que la imagen con la cual son asociados, sin 
embargo residen en un subdirectorio dierante Ovas 
extensiones son utilizadas para el Especiro de 
Frecuencias de Fourier (AMP), para las partes real (TFR) e 
imaginaria (TFI) de la Transtormada de Fourier. Varias 
extensiones son empleadas también para los archivos 
temporales. Otras extensiones empleadas son las 
siguentes: ISO para las imágenes ciasicadas con el 
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Algoritmo de Isodata, VER para la clasiicación 
supervisada, DIS para las estadisicas del Algorimo de 
Isodata y DVA para las estadísticas del ciasicador 
supervsado. 

Por medio de otra utlei, las etiquetas pueden ser 
removidas (borradas) dejando únicamente a información 
pertinente alaimagen, la que asu vez puedeser exportada 
úaotrosisiema de analisis digal de imágenes. Losmóculos 
para desplegar y manipula las imágenes se encuentran 
Contenidos en el grupo de utierias (Fg. 1). Se tenen 
disponibles cuatro modos de despliegue: MCGA con 
320:200 pireis a 256 colores simultáneos, VGA con 
404480 pixels a 16 colores simultáneos. SVGA-Paraciss 
con 640x480 pixels a 256 colores simultáneos y 
SVGA-Tseng con la misma resolución que la anterior Los 
modos SVGA son particulares a caca tarjeta de video por 
lo que esnecesano especilicar acualde elas (Paradise o 
song) se aplica. ol algorimo correspondiente de 
desplegue. No obsianie, el aigorimo que se tiene en 
SANDI os suficientemente general como para acomodar 
olras taratas SVGA una vez conocidas les direcciones de 
los puertos de cada una de ells. Varias paletas de grises 
y ono de color están prepragramadas en as Javes F2 a 
F6,lallaveF estásiempreresenvada para ayuda entinea. 
arras uerias incluyen el pertide una imagen alolargo de 
na linea o columna, el directorio de imágenes, la 
generación de una subimagen, el borrado de una imagen 
incluyendo sus etiquetas yla impresión entonos de gris de 
una imagen empleando una impresora compatdle con la 
horma HP (Henien-Packard), esta impresión es en 64 
(ivoles de gs y hasta 200x370 pixels cuda vez: para 
impúmir una imagen mayor se requiere generar. un 
mosaico, La impresora laser debe tener al menos 2MB 36 
memoria 

Para un manejo eficiente de memora, unicamente 
aquellos módulos involucrados en una transiormación 
dada son cargados en RAM. El número actual de 
programas ejecutables en SANDI es de más de 120, 
cupando todos ellos más de 25 MB de memoria, esto 
incluye transtormaciones y programas de apoyo. Las 
restricciones quo ainprevalecenen DOS cancemientesal 
ú¿reccionamientode memoria RAM. impone imiacionesal 
tamaño de los módulos, numero de archivos activos, 
tamaño de las ventanas (hasta 11x11 pixels) para: 
convoluciones, número de bandas y clases 
clasficaciones muliespectrles. Por otr lado, el uso det 
DOS permte la tansponabiidad del sstema a un gran 
"número de PCs. La aparición reciente dela versión 70 del 
TurboPascal, con el que so puede manejar el modo 
protegido del microprocesador de una PC y conel que se 
genera código de 32 bs, permi dseccionas toda la 
memoria RAM dsponibla eliminando muchas de las 
restricciones mencionadas y haciendo más rápida la 
ejecución de los programas. Actualmente se está 
¡scomplando SANDI en este nuevo marco de operación 
Finalmentatodaslas imitaciones sonun compromiso entre 
lamemoria RAM disponibl, el dreccionamiento al isco y 
alllempo de procesamiento. Por ejemplo para una imagen 
de 1024x1024 pixels, las manipulaciones al histograma 
puedontomarvarias decenas de segundos, peroparauna 
clasiicación no-superisada con sete bandas y dez 
clases, el tiempo puede ser de varas horas. Los fis en 
«el dominio espacial (por ventanas) requieren de vanos 
minutos, pero en el dominio de frecuencias el iempo 
puede ser de varias horas. Otras Iranstormaciones 





SISTEMA DEANALISIS DIGITAL DE IMAGENES, V. 20 


SELPER 


demandan generalmente de varios minutos. Estos datos 
on para una PC 36/33, 

Én cada meno se ene disponible una ayuda en Ine, 
que va de una explicación general a, instrucciones 
especficasacercadelaoperacióndeuna transformación 
Los mensajes de control y eror son manejados. por 
procedimentos especializados que residen en unidades 
de 3poyo La instalación especiica(Fig.3)es como sigue: 
odos los programes ejecutables residen en el directorio 
llamado SANDI, los manejadores para depllegue e 
impresion están en el subdirectorio DRIVERS, las 
imágenes en el subarectario IMAGES, las estadísticas de 
los datos ciesiicados en el subalractorio STATIST, los 
aschiwos de etiquetas en el subdirecorto LABELS, y los 
archivos temporales en el subdtectoro TEMPO. Aún 
cuando los archivos temporales son-automálicamente 
borrados lfndecada procedimiento que requiere suuso, 
el subarectono TEMPO es úl cuando un proceso es 
internmpido por alguna razón, en esto evento el eror 
puedesertrazadooelarchivo temporal puede serborraco 
Sin necesidad de ninguna tarea de búsqueda. El 
subdirector HELP contene todo el exo necesario para 
la ayuca en nea. la que Incluye una breve descripción de 
la operación a iniciarse e instrucciones para iniciar una 
tarea dada 











El copiado de SANDI a un disco duro a part de 
dskertes es una operación directa; los. subidiractorios 
mostrados en la igura 3 deben ser creados de antemano, 
¡en este momento no se tiene un instalador. Las imágenes 
on copiadas al subairectorio IMAGES y por medio de una 
de las utlerias, las etiquetas respectivas son creadas y 
relacionadas a cada archivo de imagen, a partr de este 
"emetass posta splcar celia ranemmación ya 
ue las etiquetas son automáticamente leidas. 


Operación de SANDI 


En el comando CASANDI> de DOS, el sistema es 
invocado tecleando simplemente SAND! esto tras al 
'desplegue la pantalla inicial de presentación y después 
¡de <ENTER> se tiene el menu principal. Tecleando el 
número asociado a la selección, o empleando las llechas 
¡elteciadoy <ENTER>, elusuano puede acceder aunado. 
las dos grandes ramas de SANDI descritas en la figura + 
Encada menú, la lavo Fl trae una subventana donde se 
da una breve explicación sobte la transformación 
invocada por el usuario, las instrucciones de operación 
también son Incluidas. Se asume no obstante, que el 
usuario tiene una lormación general en el área del 
procesamiento de imágenes y el reconocimiento de 
patrones. Estas ayudas esián disenadas para servir 
"Unicamente como guía en la operación de SANDI. 


SiLPER 


¡Siguiendo la trayectora adecuada (Fig, 1) el usuario 
aventualmente llege a la transiormación especifica 
deseada, en este punto se muestran secuencialmente 
varios submenús donde se solicita: los nombres de las. 

entrada y salida, el valor de los parámetros, 

idas poros nacer arca 
necesarias para 

lransformación está relacionada con imágenes 

0 número y nombre de las bandas son 
En aquelos módulos donde se requiero una 
¡dinámica con la imagen, se incluye un editor 
"Que permite al usuario manipular el ratón, las 






pra 
nitanamíento). La ayuda tiene también 
'SObrO el uso del editor gráfico. La última 
de cualquier meno siempre indica al inmediato 
¿Goma excepción del menú princ:palen donde se 
regresa u DOS, de esta manera el usuano puede r hacia 
“adelante o: hacia atrás en el árbol de opciones del 
dlagrama de SAND, 

limediatamento después que una transformación es 
Iniciada, aparece una Ventana donde se reporta el 

ella tarea, se proporcionan: al nombre de la 

¡alpaso y eltipo de proceso yla mea dela imagen 
bajos. Los opos ic! y Inle despogan al 
Uinal de cada Para -omposición en 
EIREIASs prericales y. el antlón Larenbo, se 
"proporcionan otros raports talas como: ol vector de mo- 
¿lala matriz de covarianza, ol porcentaje de información 
nas bandas de salida yla matriz de correlación. Para el 
rocimiento de roglones y los clasifcadores espectrales, 
ño proporcionan las estadisticas de cada clase, 
Incluyendo la distancia -M (Fichards, 1986), En cada 
CASO l situa pormio al usuario elimprimi o almacenar 
'n un arcnivo la información de estos reportes. 

El dospliegue de una imagen se hace dando su 
mbr en primera Instancia, a continuación aparece un 
menú con cuatro posiblidados de resolución: MCGA, 
VGA, SVGA Paracio y SVGA-Tsong; la solección de una 
de óstas pone la imagen en el montor. Sel tamaño de la 
magonexcedo aa do a ventana da desplegue, entonces 
o acia ol desplazamiento (crol) de ésta por mecio de 
las flechas del teciado, En la ll 
preprogramada una paleta Incal de gres, on la FS una 
Jogafímicade grises, una Arcotangente de grises entaa, 
'napaletadecolerentaF5 y para imágenes de compuesto 
.color sa eno una paleta on FG, la paleta original de la 
tareta e despliegue es lía y almacenada en a lave F7. 

'Entoda ransformación las etiquetas de la imagon de 
entrada son leidas y transmilidas al arcnivo de las 
elíquetas de la imagon de salida, el contador sobre el 
número de procesamientos es Incramentado en uno y la 
liqueta respectiva de historia de procesamiento es 
añadida, Todas las operaciones de convolución están 
rostringidas a ventanas cuadradas, impares de tamaños 
“que van desde 3x3 hasta 11x11 pits. Alrededor de 800. 
KB de memoria RAM es requerida para correr cualquiera 
“plcación de SANDI, pero s por alguna razón la memoria 
ds escasa, entonces todos los casifcadores pueden ser 
nvocadosseparaciamente ya que éstos son autónomosen 
cuanto al despliegue y manejo de la inagen. 
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Aplicaciones 


SANO su sena de ani gu! e iagures 
de propósito general aplicable a cualquier 
Investor datos seria a rocio da 
imágenes y el reconocimiento de patrones, len hasta 
Srorasiéntasson uusohascoentapercepcióntemoa 
También SANDI es empleado para apoyar las labores de 
enseñanza en un curso sobro procesamiento igilal de 
imágenes que se ofrece regularmente a nivel de maestría. 
El sistema es empleado en prácticas de laboratorio para 
lurtar a los estuclates el potencial de los métodos 
expuestos en el cuso. En acividades de Investigación 
SAND! ha sido empleado con éxito en proyectos para 
hacer análisis de texturas en imágenes de radar de 
apertura snética, en clasficaciones espectrales de 
imágenes multespectrales de satólite para dilucidar uso 
de suelo, en segmentación de imágenes multiespectrales. 
[cir enero abra ei 
urbanas y suburbanas, en estudios 

cuencas cerradas y en estudios morfológicos de 
partículas ullizando Imagenes de microscopio 
elctronco 











Conclusiones 


Homos giseñado, desarrollado o implantado un 
paquete para andliws digtal de imágenes, de caracter 
“amigable y operado por monús. El paquete funciona bajo 
DOS que os un sistema operativo de ampli uso poro quo 
"o ransportablidad está asegurada en una gran cantidad 
do microcomputadoras. Con a aparición del TurboPascal 
7. con el cual es posible gonorar códigos de 32-bis y 
¿Sreccionar oda la memoria RAN dispontbo,os médulos 
e SAND! operarán en forma más afcionte y so podrán 
mavejat imágenes de hasta 204842048. Por oro lado us 
importante puntualizar que con la tuación actual do ls 
Imicrocomputadoras el andisis de más de una Imagen 
"mullespectral de satóito es poco realista en esto 1po de 
plataformas: una estación de trabajo sería adecuada para 
' procesamiento de varias do estas imágenes, El plantear 
un sistema transportable impone nocosariamonto 
limitaciones en cuanto ala capacidad de análisis: requerir 
platalormas mas poderosas harian monos transportable ol 
sistema, En este contexto, el alcance y tipo de proyectos 
que pueden atacarse con un sistema como el presentado 
aquí deben de ser cuidadosamente definidos de 
antemano. Proyecios de pequena y mediana oscala son 
factibles con SANDI, los proyectos de gran envergadura 
deben ser trabajados con estaciones de trabajo o con 
¡grandes computadoras. No obstante, es importas tener 
en mente que el desamollo delas microcomputadoras os 
notable y que en un futuro próximo su capacidad será 
comparable a la de una estación de trabajo, de esta 
manera la evolución de SANDI a nuevos complladores y 
«istomasoperaiwos para estas plataformas es básico para 
mantener e incrementar su potencial de análisis. 
Finalmente, deseamos expresar. que SANO se 
encuentra disponible, Ire de costo, a la comunidad 
especialistas dedicados a la percepción remota y el 
análisis digital de imágenes Para mayores detalles 
¡Comunicarse con el primer autor (4 L) 
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Abstract 


Version 1.0011he Digital Image Analysis System (SANDI), previously reported (Lira.tal,, 
1989) has boen redesigned and reedited. Actual version. 2.0 carry more than 70 modules. 
úoperated by menus. All modules are written in TurboPascal 6.0 operating under MS DOS 5.0. 
One MB RAM and mathematical coprocessor are required besides a hard disk to store tho 
digital images. The systom, working on a PC-386 or 486 ¡s divided in to main groups: 
Enhancements and Processings, and Pattern Recognition. A group ot utiles display and 
'manipulate images ¡s available in several locations ot the package. The display window is 
edito scrollby preprogrammed keys with several colour palletos including black and who. 
Images are displayed through commercial video boards and are stored in free MS DOS 
format related to control labels that includo information about size imago, history of 
processing and a comment by the user, Main applications are devoted to Remote Sensing. 





INTRODUCTION 


igh productivity rato in now methods tor 
digital mago precossng and paro recog 
lion has been sustained in the past decade 
Actual ronda show no sign ol decis, cn tha 

Jalizod journals and astcias see tna 
lo pace, Hence, appears tna abroad 
(adosrablato1ostnowy designed modetator dgtalimago 
analysis. Tho environment should traduces into general 
purpose mocilar sysiomwnich bos portables possible 
Tha great development in tho microcomputer held. com- 
nad with now gh resolution commercial video boards, 
suggests that a PC-386 or 486 may bo an adoquato 
platirm lor such a system. The system should be easy 
Vpgradable since s main objectve sto support research 
cies inigtal Image analysis. Eacn ndcion or updat 
ingota module increases the power analysis! Ine system, 
bringing fngobanatts to te process as providing capa: 
lies orfutnary and commercial applications. Imposing 
“APC based system restricts he use ol cen operating 
systems and delnos a number l languages to bo em. 
ployed, However, he progress in spee, memory devel- 
pmentand operating system for PC. predicts o avai- 
abliy an ample environment Tor image processing in a 
lersocablo future. 

"Avaraty of commercial packages is avañable (PC, 
IPLIB, MIPS, ÍTEX, VISILOG. IMAGE-PRO. ERDAS, 
DRAGON)" ranging inpricestromone thousand to several 
tens ol ióusande ol dollars, and wortang under various 
latorms ana environments. Some olher packages are 
Publishad ín the Kterature win varied capabsites and 
environments ol operatcn (Bemstin and Fermesmough, 











Tri lo al 1 ar an endorsment 1 rose brands 


1978; Wharton and Tumer, 1081; Landy ot al, 1984 
Harmngion et al. 1986; Eiyant, 1990; Sader and Wine, 
1991) including algonthme for graphics and image pro: 
cossing (Pavidis, 1982). Soma ol them aro spocializad 
lowards a specit task (Wnaron and Tumer, 1981: Bryant, 
1990), andothersare ortertodtowarde a specili held such 
as Remote Sensing (Bernstoin and Fermoyhough, 1975; 
Harmngion et al. 1986), an expert aystom (Schowengerdl 
and Wang, 1989), or for education (Sader and Winno, 
1991), 

Vio present hare, a general purpose modular system 
operated by means ol menus, portable under MS DOS 5.0, 
'assiyexpandable and fuly wrten in TurboPascal 6.0 us 
intended to support actils inresearch and davelopment 
in the area of digital image analysis, evenbough lt is 
'omployed intecting as well The structure and operation ol 
this system is described in the folowing Iwo sections, 





STRUCTURE OF SANDI 


Before quanting the propertes ol-a pattern 1 is 
usual required to enhance Ít. A pattern or a class. ol 
patterns are notin general well ierentiatec with respect 
to the rest of the imago, Noiso, blurring, intorpixa! rola- 
ú"onshios or lumination innomogensitias may obscura or 
hiden the presence ol a pattern, rendering dificult is 
recognition. An entancement or processing is in general 
requieg lo put n evidence tha paltem previous lo any 
recognition procedure. This basic assumption crives thó 
main structure (Fig. 1) of SANDI. Processings and En- 
úhancements are classiled according lo the type ot op: 
raton they employ, thus, Radiomatric Enhancements 
stands for manipulatioos to the histogram. Geometric En- 
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hancemants aa to put ín avidence some gesmevie as» 
pects of the imago. Fidelty Ennancements aro basically 
festorations to the image since they tend to approach, as 
much as possible, a degraded image to an ideal modal. 
Feconsioion and Compression tha! utiizes zoom by 
pines, and reducion of tho Image by elíminating one lo 
ch hor pix. And Ensemble Operations to transtorm 
ha image as vola, Each ol those sections Include a etot 
operations according o tha diagram of SAND! (Fig. 1) 
The secton ol Pattern Recogaan is subdivided nto 
Spectra and Spatial patterns, eventiough those sections 
are very much rolatod, Segmentation ol spactral patterns 
include technique ol ens slcing employing the ni 
logram, a non-supervisod classifar ulizing tha Iso 
algorihm (Tou and (Gonzalez, 1974) and a: bayesian 
¿lasolier fichards, 1986) using a maximum Iklligod 
alo, A togion growing algoritm (Wah, 1987) a included 
wit tha Tolowing uniormiy prodicatos: average car 
ne, corelalzn and antropy based an co-occurenco 
matrix. A contextual classlier based on co-cecurrenco 
matrix and entropy orlterion works lor monospectral im+ 
gos and 1 special sutod for radar imagos analyais, On 
he side of spatial patterns, there is segmentation oriented: 
towards lines detection, algobraic operations among 
bands, pet operations tor clansifiod images, and math» 
emalical morphology (Serra, 1982) on 2D and 3D with. 
selectod applications. A set of utilities 10 support the 
manpulatons o! magos la Included n selected places ol 
tha system 
Intermaly SANO! a distributoc in programs and unit; 
lis menús form a set apart ol programe. Units group a 
«number ol procedures tor specialized applications, some. 
¿Xherunis ncludo procedures or lies, menus, reports and 
display handling. Addíition of new algorithims involves only 
the updatig ol a menu and the uni whore tho program 
(procedure) ls to be included, recompllation Iherefore is 
costicied 10 is part ol no syatem. Conwol ol manus 
roparte on erors, and about ho progress al an operation 
fs carmod out by spocial purpose dosigned procoduras 
lhat can ously bo modiod o accomodato now tasks 
Iímagos aro siored In tres DOS format win block aízas 
“quel to the number of columns in the image, labels ara 
hancled in soparato los ol 128 byios blocks. Maximum 
eres to 1024 1024, basically dueto 
rmiations. Fourier Tranaloram in med 1o 2 order 
nd region growing to 512 x 512 sizo. Spectral 
classification, both supervised and non-supervised 
peralo lor 7 bands ana 10 classos maximum. Canonical 
analyals n he olhethand works for up lo 5 classes and 7 
Bande, A as malfxranaposo algoritm (Huang, 1981) a 
amployad tor FFT when the image is 1024x1024 size, 
'One byte per pixel ls employed in every image with the 
excoptonol temporal nos whoro varying degroo af prect- 
son May require Iwo, four or elght bytes per pixel Tem- 
poral fos ara automatica erased once an operation ia 
ended, Control labels are annexed to label files (Fig. 2) 
employing one of the utilities. Labels include the size ol the 
Image; pixels per line and, number of Ines, number ot 
processings received, a comment by Ihe user of 80 char- 
Aclors maximum anda comment on the history of process: 
ingisatided each time an operations applied to ne image 
Where k is an Integer, On: the other hand Im = numbar of 
history labels, then lfm (mod 2)=0, we have Ihat k = m (div. 
2)+ 3, lhenvso,k=m (dw 2 ).+2. 
"magos aro idonted by thor name aná no oxtension 
'POI, those with no labels carry the extension IMA. Label 
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ios havorh samo name and extension trino imago they 
are associated to, but ihoy rosido: ina diforonto 
«ubdlrectory Other extensions ar uso or th Irequency 
Fouiór Specirum (AMP), real (TER) an imaginary (FI) 
parts o Fourier Transtorm. Various extensions aro used or 
temporal fos as wok. 

Other oxtonsiona aro ¡ho folowino; 1SO tor 
or Supervised ciasiliaa data, VER or supervise 
silecdata, DIS or]sodata Stataucs, and OVA for Bayesian 
Staisics. By means ol anoiher uy, abeto may bo 
¡moved (cralhod) leaving on te Information portaining 
10100 Ímago, wieh ín tu may be oxported 19 somo othor 
ystemol analysis The mogulos to display and manipulata 
imvagos aro contained in tho Uites set. Theo display 
modos aro avallblo. MEGA, wih 320x200 plsols at 256 
Colors, VGA wi 640x480 plxls a 16 colors and SVGA 
wi GAOK48O al 256 colors is Inst mo a only val lor 
Paraiso video board, however ho algorihm is gonaral 
'ónough lo accomodate other boards, Sevoral paltas, n 
gray and color tones ara proprogrammed in th koya F2 
FS, o F1 koy la roservod for orina holp, Other villas 
incluco tha prof ot an imag along any no or column, 
diroctorataof magos, anda gray tonos prato!an mago to 
' lasarjot printer compatible wih HP standard norm, 

For an oliciont managomont of memory, any those. 
modules involved in a given translormañon aro load nto 
RAM, Theaciual number l oxecutablo programs in SANO! 
ls over 120 occupyng more than 2 MB ol memory, this 
includes transtormaions and support programs. Resto 
tons til prevaling on DOS concomingIha adetessing al 
RAM impose límiations on modue izas, number les, 
sizool windows forconvolutions,andnumber ofbaands.and 
¿lasses In mutispectrl images and speciral iasaifica- 
tons. On ho other han, use o/ DOS permi he porabity 
ofihosyalemtoalargonumberolPCS, Rocentappesrence 
o Pascal compler compatible with TuboPascal and 
running, under OS/2 may be one ol the posibiites ol 
development for ino naar futuro, anothor atematwa ls 1o 
'employ TurboPascal or Windows, This Imation a 
fact a trado ollo avoid oxcosive dk addressing and to 
ú'nssuo a reasonablo processing time. For example forn 
imago 1024x1024 size histagram manipulaions may tae 
savoraltens of seconds, bu non-supervised clasaficaion 
vit seven bands and len classes may lako several hours. 
Flors ín spatial domaln require several minutes, ut in 
frequency comain up to several hours. Other translorma: 
tons usually spend sevoral minutes, Thoso figures are for 
2386/83 PO. 

An integrated tutorial is accesible rom esch meno, 
¡golng from general axplanations 19 specilisinsructons 
aboulpartculer operation ol given transformaton, Con: 
trolmessagesancerrorreporsarohandied by specialized. 
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procedures resicing in support units The specili installa- 
on (Fig. 3) ls as follows al execulable programs reside in 


A director named SAND), vers to display and prin 
Se in 6 aubalociory 'DAIVERS, Images are ln e 


subdirectoty IMAGES, starstos for classic data are in 
thosubdiectory STATIST, label les are intho subdirectory 
LABELS, andiemporalfibs are in the subdiractory TEMPO, 
Evenihough temporal Nes are automatically erase ar the 
end ol overy procedure that requita Ineir uso, tha 
subairectory TEMPO, is useful when a process is Inter- 
rupted lor any reason, ln is event, an error may ba traced 
lup ora lemporál le may be scratched wiholt any search 
task, The subdirectory HELP contains all ho text matenal 
lor orlne help which Inclucde a. Brie! deseríption of the 
úoperatiónto'bo undertakan and instructions to start givon 
10sK: 













SANO! 
— DRIVERS: 

|— IMAGES 

|— sransr 
[|— LageLs 
|— TEMPO 
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Copy ol. SANDI from diskattes 10 a hard: disk is a 
ústraightlorward operation, Ine subdirectoris shown in ig- 
ro 3 must be created in advance. Ar his timo Meca ls no 
program for autormati installation, Images are copied to 
ho IMAGES subiioctory, and by means of one ot the: 
Vblitios, labels are croatoc and annexad (relator) to vach 
image Ilo, at his mo any transformation may appled: 
inco labals aro automaticaly road 


OPERATION OF SAND 


AN iheprampt CASANDI>0l DOS tro systom¡sinvoked 
by1ypIngSANDI,ihis brings ho main menuinto display, By 
using Ito mouse, the number al selection; or the-arrows ol 
ha koyboard and <Rotum», tha user may accoss ono ol 
thotwobranchasdepicted nfIgure 1. Atoach menu he key 
Ft brings a subindowwhero abrielosplanation abone 
transformations Invoked by ho usar ls gisplayod; instruc+ 
tons for each translormation aro included It is assumed 
however Ina! Ino user has a genoral background on digital 
Image processing and pattern recogritlon. This tutorial 
Intonded as a guldance to the operation ol Ihe system. 

By following the right pata (Fig. 1) the usar aventually 
arrives at an specifi transtormation, at is point several 
úsúbmenus are sequencialy shown; the names o! the input 
and output images, he valve ol parameters. nvolod 
whoserangeol valdiy is writen between parenthasis, and 
somecther decisionstostart tha application, the transtor- 
matlon¡srelated 1o multispectral images, the number and 
name. the bands are requested. Inthose modules whore 
“dynamic interaction wlh Ine image is required, a graphic 
cir ls incluced that enables Iho usar to manipulate Iho 
mouse, tha arrow keys and palettes. This ls useful lor 
classilersandcanonicalanalysisthat neecttho definition! 
training ido, The tutorial inclucies detallsaboutthe uso ol 
the graphic editoras well. haastoption ina menu always 
ndicates ho previous one, with the exception of the main 
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menu wbich toturys the user to DOS, n ni way tha usar 
vay 90 forward or backward in tha tras o options 

"As soon as atranstormation stated, a windowwth a 
reporton he progress le ask appegrs ho name o! he: 
mago, (ho slep andl type of process, and ina Ine under 
analysis ara given. The nia and final timo aro cisplayed. 
altho oncol every task. For principal components iecom- 
úpositonand canonicalamnalysisolharreporis are provided 
such as mean vector, covarianco matrix, percomago al 
information in output bands and corelaion matrix. For 
region growing and spectraclassiers class statistics ro 
providod, nclucing ho J«M aistanco (Richards. 1986) lor 
supervied classification In every case th system allows 
tiousorta printonhardcopyarstore nal maintormation 
finos reports 

“he display ol an image 1 dono by typing ts name, 
altorwardo, a manu with vos rosolulione ls brought 
MCGA, VGA, SVGA; tha. selection ot ono ol those 
immediely put the imago onto tho monitor screen, ll tho 
Image sizo oxcoeda the display window, serall a then 
“avalale by means ol he arrow keys ol o keyboard. A 
black and whit palete 1s prestored in tha: F2 hay, 
Logaritmio gray tonos in F3, an Arctangent n A, color 
paleto ln FS, ana for color composta ima 
Ín FS; ho original palena,ol ino video bogr 
storod in F7 

In ovory transformation, mago | 
transmicion oupuA label lato mago, tna counter 
¡on tho number ol procossings is increased by one, anditho 
ho rospactiva history of processing labo! Is. acidos AN 
coNVOlulonoperatonsaarorestrcieita cd windowsrang 
Ing insizosfrom3x31o 11%, About S00KBO! (roo RAMaro 
hoaded to run any application from SANDI, bu or somo. 
roason RAM memory ls low, al classllers may bo run 
Soparatoly from DOS, Tho reason for is sta clasara 
aro somewhat autonomous including display capable 





















APLICATIONS 


SAND! ls a Guonoral purpose digital mago analysis 
systom applicable to any aro ol rosiarch where pattern 
recognition nd image processing are required WN em- 
phasison Remote Sensing applications. Actually, SANDIIs 
úlizad to support a course on Ihese opics that we regu: 
1ariy offer at mastor degree level In Ine area ol computing, 
sciences, The systom ls Usod as a basis for laboratory 
practicos in a hands on experience lor students, In ra- 
search activiles, SANDIhas boon succestuly employed in 
projects for textura analysis of SAR images, ln 
mullespaciralclassificalon o! satelite images forland uso. 
purposes, In segmentalion ol muliespectral images for 
Iinologies identiicaton, in segmentanion of urban and 
úsuburban areas, in geohicrologíc studies ol closed basins, 
nd in morphologíe studies of partiles unlizing olectronic 
microscope images. 











CONCLUSIONS 


Wo havo designed, developed and implemented a 
user Inendi and menu ven package for digital maga 
analysis. This system operates on a widely implemented 
operating system as DOS is, however, a more powerful 
environment ls required to Increase the capabilty ol Ihe 
package; OS/2 might be one altemativo, For practical 
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purposes mulliespectral. saven bans images of 
1024x1024 size represent Ine real lim for gia) image: 
analysis on actual PCs. Biger sizes or more bands would. 
requíe a Workstation, Since the idea is lo havo.a portable 
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present here. Small and madium scale, ana research 
projects are certaint feasible wit is 001 The analysis l, 
let us say a whole framo¡ot a LANDSAT TM image is not 
reasonable by a system runping on a PC Ín he actual state 


liber or ese machines; the near future looks however 
very promising. 


system willy minimum hardware requirements, a compro- 
miss must be made concerning the volume ol data to be: 
angled. Thíshas implicallons orina ype and scopeoftne. 
Projectsthat can bo worked out by a systom lko tna ono we 
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TELEDETECCION DE LA ERUPCION 
DEL 19 DE ABRIL DE 1991 EN 
LA ISLA FERNANDINA 





Hugo Pérez Oviedo 
CURSEN. Ap. 8216: Quito Ecuador 


Introducción 


El presente articulo describe la información vulcanológica obtenida en la Isla 
Fernandina, utilizando datos de sensores remotos, principalmente del satélite LANDSAT 
(Sensor Mapeador Temático), que fuera recibida, grabada y procesada por la Estación 
Cotopaxi del CLIRSEN en diferentes fechas. 

"Observaciones de campo, efectuadas en septiembre de 1989, y un estudio realizado.con 
varias imágenes LANDSAT, una fotografía de cámara de formato largo y una imagen SPOT, 
permiten aumentar el conocimiento dela evolución volcánica reciente de a sia Fernandina. 

Trabajos previos, con datos de sensores remotos sobre esta isla, permitieron obtener 
mapas topográficos y cartografiar flujos de lava, fisuras y conos (Munro y Mouginis-Mark, 
1990). Una separación de los flujos por su edad relativa no se pudo conseguir por falta de 
contraste espectral (Munro et l,, 1991). 

Este volcán, La Cumbre, es el más activo del Archipiélago de Colón, con 21 erupciones. 
registradas desde 1813 (Simkin et al, 1981), 11 de las cuales corresponden alos últimos 33 
“años. Luego del colapso de caldera en mayo de 1968 (Simkin y Howard, 1970), las siete 
últimas erupciones se han efectuado a través de fracturas circunferenciales íntimamente 
relacionadas con la caldera. Con excepción de una, las erupciones son intra-celdera y, 
aparentemente, pertenecen a un período de relleno de la misma. 

Este proceso de relleno parece ser repetitivo, como lo sugieren las intercalaciones 
úhorizontales de lavas y tetras que se ven en las paredes de las mesas NNO y SE; expuestas 
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luego del colapso de 1968. 





EVOLUCION VOLCANICA ACTUAL 





11 1972 y 1977 fueron erupcionados dos jos 
¿de lava sobre la mesa SE dela caldera, y en 
1973 ocurrió uno muy corto en a pared ESE 
vecina ala misma mesa (Fg. 1) 

El tujo de lava de agosto de 1978 tua expulsado por 
tna fisura ciscunterencial entre a pared yla mesa de! norte 
de la caldera, cubriéndola parcialmente para descender 
en cascada hasta a base de la caldera. El uo de lava de 
1982 saló por varios ventos relacionados con dos falas 
etrcunterenciales ubicadas en la planicia de la cima, en el 
borde sur de la caldera (Mouginis-Mark et al, 1990). La 
lavaesdeltipo pañoehoe que varia itipoa-2al descender 
Una corta pendiente. El fjo de lava intra-caldera de 
septiembre de 1988 estuvo acompañado por una corta 
úoxpulsión de escoria, la misma que alcanzo la planicie de 
la cima (zona B en Fig_1) Este lujo salió por una fisura 
localizada junto a la pared suroriental, descendiendo 
hasta la base de la caldera, que en ese entonces estaba 
parcialmente ocupada por un lago, que quedo reducido a 
lina pequeña laguna luego de la erupción 

Observaciones de campo sobre la planicie deta cima, 
revelaron la presencia devarías tallas circunterenciales de 
coro desplazamiento normal, las mismas que son 
aprovechadas para las erupciones más jóvenes, como la 


























de 1982. Estas tálas están inimamente relacionadas con 
la evolución de la caldera y podrian explicar la morfología 
plana tipica de las cimas de los volcanes do Galápagos, 
Éstas falas también pueden hospedar digues como ingk- 
¡can Me Bimey y Véllams (1969) y Nodlle (1973) en sus 
explicaciones sobre la morfología de fas Islas Galápagos, 
Continuos desizamientos en la pared oriental de la. 
caldera provocaron una migración de su borde desde 
1984 (zona C en Fig. 1). Esto se pudo observar 
¡comparando una lotograta de cámara de formato largo, 
acquitda enoctubre de 1984, y la magen SPOT del 28 de 
Junio de 1988. La zona B corresponde al colapso ocurrido 
"durante la erupción de 1988 (Boletín SEAN, 1988). 


ERUPCION DE ABRIL DE 1991 


Después de ocurrida la erupción del 19 de abril, sólo 
816. de mayo de 1991 se pudo obtener una imagen de la 
caldera de la Isla Femandina sin nubes (Foto 1) 

Una composición a color de esta imagen, con las. 
bandas 4, 5 y 7 (roo, verde y azul respectivamente), 
"permitió observar dos aspectos importantes de la úlima 
erupción. 


Dos nuevos Mujos de lava hablan sido erupcionados 
¿entro de la caldera. El principal habla salido por el 
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ventoda 1982, atravós deunaisuracircunterencial de. 

la pared oriental, 950 metros al norto ds la mesa ESE. 

on la, zona donde en 1989 se notaron continuos. 

derrumbes, Sus coordenadas son 91* 3040" LWy 0% 

29'24* LS. Luego descendi formando un canal de 

2,18. km. de longhud hasta la baso de la caldera, 

donde cubrió parcialmente el flujo de septiembre de: 

1988, Deucuerdoaldreavcupada oral fujose puede 

dedudir que es de mayor volumen que los 

anteriormente cilados (Fig. 1). El Mujo más paqueño 
había salido por una fisura circunferencia! del nora, 
por dondo ya escurió una orupción en 1978, su 

Tonglud es ce 850 mts. 

+ Haciendo una comparación de la imagen del 18 de 
mayo de 1991 con las Ímagenes del 5 de enero de 
1090 y del 13 de marzo de 1991 (Foto 2) tocas coma 
composición a color citada, se puedo observar 
elaramente un cambio.enatonalidady enta textura de. 
ll caldera (excepto en el área ocupada por el nuevo. 
JO). Ena imagen del 16 de mayo, latonalidaces más. 
élara y la toxtura es más suave, indicando que la 
caldera fue cubiena por un marto de cenizas, Esta. 
Misma observación se tiene de la paria noroccidental 
de la cima y en los flancos occidentales del volcán, 
hacia donde se extendió la nube de cenzas. De 
acuerdo a reportes de: dos testigos (Bolerín. GVN, 
1991), la expulsión de cenizas acompaño 
Intormitontemente a casi toda la erupción. notándose 
“que en ol Cabo Hammond se acumularon 5 - 10 mem. 
de cenizas. 

Comparando las mismas imágenes, se puede notar 
nmbién unoscurecimiento ds una pequeña zona vocina al 
¡fuevo venta (A en Fig. 1). S/ comparamos estas imagenes. 
¡con las del 20 dle marzo de 1985 o con las del 28 de junio 
de 1988 del satélite SPOT. se nota el mismo 
'oscurecimiento en la zona cercana al sur del vento de 
septiembre de 1988 (0 en Fi. 1). Esto.se interpretó como. 











¡una acumulación de escoria, contimándose luego con el 
Boletín GVN (1991) que reporia la acumulación de un em. 
in esta zona, 


ESTIMACION DE TEMPERATURA 


Al roakrar la combinación color de as bandas 4, 5y 
7 con colores an. verde y ojo tespecvamente, se pudo 
ota que varos pelos iocalzados muy corcadel vento y 
ne(canaldelava  aprorimacamente 700 mis. el vento, 
parecieron con coo rojo, Dentro dell, otos piselos 
poseen los valores más altos de intensidad (37) ena franja 
Errar, a misma que dere una resolución do 30 mis. 

Usando las técnicas tormalos de esimación de 
temporaturas apart do datos atoltaris, quese bosanon 
lla apicación de la Ley de Radiación de Planck (Francis y 
Rotary 1987; Rory eta, 1988; Glazaotal 1980).y con 
las tablas de equwalencia eme números digilos y 
lemperaturas para las bandas 4, 5, 6 y 7 de LANDSAT 5 
elaboradas en el CLIRSEN (Pérez 1990), se establece una 
estmación de 196'C para estos pieles. Es dect, que 22 
¿as de terminada a erupción, exsttan todavia nas dal 
ll, de por lo menos SO x SO mis, con erparalras 
promedo de casi 200€ 





CONCLUSIONES 


La información satelana, recibida y procesada por la 
Estación Cotopaxi del CURSEN, ha permilido adult uo 
conocimiento rápido y casi completo de la erupción del 
volcán La Cumbre de la isia Femnandina, ocurrida entre el 
19 y el 24 de at de 1991 

Esta erupción se acompañó de fases explosivas que 
úyectaron una buena cantidad de cenizas, las mismas que 
cubrieron la caldera. la cima y los flancos occidentalos. 
También expulsó escora que alcanzo el borde. de la 
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deraentacinajunto alvento, Loslujos delava salieron — AGRADECIMIENTOS. 
por fisuras circunterenciales de la caldera. llegando el 
mayor a cubrir casi todo el londo de la misma, 
días después de terminada la erupción 

nian partes calientes del lujo (cas 200*C) 
detectables con el satélite LANOSAT en el rango in 
del espectro electromagnético, 

Desde 1984 hastalafecha, senotaun cambio drástico 
enla lorma dela caldéra por os desizamientos y colan 

alzados. 





Deseo dejar constancia de mi agradecimiento a tas 
dades del CURSEN, por autorizar la realización y 
ación del presente trabajo, y también l personal de 

cesamiento y Producción, en especial al Ing, Lino 

/erduga, por su Colaboración, prestada durante el 

umiento de las imagenes, 
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FOTO 1: IMAGEN DE LA CALDERA DEL VOLCAN 
LA CUMBRE - ISLA FERNANDINA. 
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Abstract 

'Synthetic Aperture Radar (SAR) images are in general monospectral providing a rich 
content of scene texture. Due to the conerent electromagnetic radiation employed in 
ú¡generating SAR images a grain-lke speckle noise is superimposed on them that obscures 
ño toxture information. 

In this research, a method to detect (segment) texture objects employing a region 
growing schemo with specific niformitypredicates designed for this ask is presented, For 
An overall texture segmentation ol the whole image, a hybrid approach is proposed. An 
entropy eriterion, based on the co-occurrence matrix, is applied, creating a synthetic band: 
which ls Incorporated with the original, into a two band multispectral set. A clustering 
analysis with an Isodata Algorithm is done over this set yielding a map of textures. An aerial 
and a satellt SAR images are used to exemply the methodology developed in this work. 

A discussion about speckle reduction is given as well providing details about its impact 
on texture segmentation. Speckle filtering is applied to both images before texture analysis. 


ls worked out. 





1. INTRODUCTION 


ynihotis Apertura Radar (SAR) imagos prosent 
'nawadays a great potential ot applications in 
manyareasol research (Kolawolo, 1991: Kwok 
tal, 1991; Marom and Shemer, 1991; Leber! 
1990, Wang and Ho, 1990). The information contento! SAR! 
images is complementary to that provided by optical sen- 
sors such as LANDSAT or SPOT. Being an actwe system, 
SAR is ablo to uminate tho target area under iftorentook 
ingles, Various polarization schemes and selected wave- 
lengihs. The. ampltude ol Ine backscattered raciation 
depends basically upon the roughness and siope of the 
target, he incidence angle, the wavelength, and dielectric 
volumetfic constant which in tum is a función ol the 
Physical properties o! ihotarger. SAR images show a great 
detal! about the target surface, Obscured ont by an effect 
known as specke 
SAR systems generato a wavetrainot coherent phased 
úlectramagniaio pulses inatinteract with ne surface under 
study. The geometric characteristics of tha wavetraín 
(Elachi, 1988) are sucn that, a resolution or unk cel ls 
iluminated by the radiation and ne intensity of the 
backscatored pulses determine Ihe pixel value assoct- 
ted 9 a cell However, in he coliIhere are a number ol 
úobjacis wilh random shapes and orientatons Ihat act as 
ú¡calterers when the size ¡sol the order of he wavelength 
úfincidentraciation, Those objects backscalteriheincom- 
Ing wavetralnin diferent directions and phases. As aresull 
úoLihis, an interlerence phenamen ocous at generales. 
constructivo and destructiva interaction ol radiaton, The 





interforence phenomenon is known as spockle conerence 
lloct,simsiar in apporeance in tha image as ino graln ota 
Photographic emulsion. Asa definition, specido (Dewaale, 
'Wambacg and Costoriinck, 1989) stho resulto construc: 
ve or destructivo Interterence ol coneront wavolats 
dephased between each other, backscatiored from a 
surface whose roughness is comparable to he wavelenght 
ot employed radiation. As such, speckle is dependent on 
tte local texture ol ha scene being studied. n this paper 
a seres ol results to segment textures in SAR images, 
'employig a region growing scherme for detectión ol spa: 
if texturo objects, and an entropy measurement, based 
úoniheco-occurrence matrix, lo generate amap oltextures 
¡recty related to scene surface aro presented. 


2. REDUCTION OF SPECKLE NOISE 
IN SAR IMAGES 


Oneway to physical reduce the spockle and random 
úholse 5 by inconerent average ol N-looks (N = 4 for 
SEASAT) ot the same scene; is reduces however tho 
“azimuth resoluonol he image as well Inconerentavorag- 
ing s the computation ol tr arithrmotic mean of N-looks, 
this kind ol average smear the image and reduces the 
údetalis and requires several observations ol ho scene 
Which aso not always available, 

Thecther altenatveo! speckle suppresionis y digital 
fiteriog the image, this demands some assumptions about 
Specs behaviour, ForoneJook intensity images, Ine cien- 
si function ol the speckde fed s described by a negative 
úexponental funcion, As N-look intensty images aro con- 
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sidored (Leberl. 1990) the density función is ratner a 
¡Chi-Squared with 2N degrees of freedom, as N increases. 
Ihis function approaches a normal cistribution. Theoreti- 
caly (Elachi, 1988), specide noise tor one-icok images 
'obeysa negative exponentíl and hence is multpicatve in 
he sense Mhal ts standard deviatica is proportional1o ls 
mean (Lee, 1989) However, has been argued (Frost and 
Yurovsky, 1985; Tur et al, 1982) that 1ne muliplicatve 
models not ahways correct. On the other hand, species 
ot always considered as luly developed (Nezry, Lopes 
¡and Touz), 1991), since lor urban areas in high resolution 
magos, low scalterera are included ín he unit cal 

In particular, SEASAT and ERS-1 images are Nulook 
averago, and experimental rasuls (Lee, 1981) show that 
te density function ls. a. gaussian Experimental results 
(Dewaelo, Wambac and Oosterinck, 1989) show as well 
thatín SAR ímages, Ino pixel have a sight correlation wit 
their closests neighbors and practicaly faces away. lor 
second or rd neíghbors; as a result of this, he specide 
noise may be assumed lo, be while. Eventhough tha 
úspeckle ls theorelicaly very mucn corretaed, in practice 
ho near whit noisa assumption for SEASAT. or ERS-1 
Images may be the resulto! Nuock average. line corrola- 
ton becomes naglígible a local mean and variance algo- 
film corrocióón must be dorivod, and since apeckle la 
directly related lo local. texture: the algorim must. be: 
adapive. 

For a non-correlated space with negativa exponen- 
al density function a Correction model was developed 
(Kuan, 1987) that load to a so called maximum a postenor 
(MAP) ter. Another model, based on experimental data 
(U, 1988), follows a non-white gauasian distribution rea- 
¡Sonably well o Speckle histograms. Based on correlation 
results betweon nelghboring pixels, a signal dependent 
correlated specklo model ís adopted. The speckle sizo 
¿corrolaion, ls in is model, approximately he size of ho 
Incoherant point spread function e the system. 

Tho spacklo models make, in one way or anctner, 
several simplíicallons (Lober. 1990: Dexragie, Wambacg 
nd Oosterlínck, 1980) that in practica do not work: 

1 The incoming alectromagnetic waves pertecty polar- 
izo 

2- Tho Incoming lectromagnetio wave is pertactlymono- 
chromatc. 

3- Tho number ol scattorars in a unit coll is necesarily 
large. 

4.- The number of scatterers in a unit cell ls constant. 
instead ts a random variablo, 

5- Tho phases associated to Me scaterors do have a 

¡uniform density function, 

"Tnís means that speckde model based fiters do not 
portorm quite well as Inoy dovate rom those facts. In 
cicition to tw, ter lend to blurr the edges of the image 
since they always introduce some kind oflow pass fiterng. 
“The goalistoreduce the noise. preserving in most degree 
ho ediges, Ihe borders and ns in he image: The type ol 
fiter foc Iris. depends uson the density function or Ihe 
¡geometric structure al the speckle; mora detal tri are 
given in the next section. 


3. DIGITAL FILTERING OF SPECKLE 
Begining vih a: simple approach, several acapive: 


¡lters were applled to two selected (Figs. 1 and 2) SAR 
mages: bolhaerialarrisatelite, Assumingthatintomoge- 
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¡nous regions the pieeis population presents a unilormly 
randomvalves distribution, and thatin edge areas the pixel 
population belongs most lo a restricted and welghed. 
average over a window was applied, bis was the casa lar 
he kenearest nelghbor and sigma fiters. Having applied 
nose files oc endows ol $ 5. point Ike ncise remains, 
so a cascade wit a median fiter seems appropriate. ne 
"cascade combination Sigma + Mein, for vncows 5 x 5 
ano 3 x 3 respectively. Yields reasonabl good results, 
ú'venthough the edges are some! blured, 

Ina second step and assuming tal that speckle ls an 
uncorrelateg: gaussian distributed density function and 
mutipicalwe 1ype noise leads to the application o Les! 
mutipicatwe fer (Leo, 1981) wih a 5 x 5 window and a 
noise variance 0 7 = 026, whose valuo ís lakon from 
'expermental ata, e. Irom ine imago sell. Results o vs 
ter show a sight visual improvement wi respeto ho 
alorementioned cascado combination. WN a signal der 
pendent correlated speckle model a generalzaton ol 
Los'smutiptcativafieris produced Inti generalzaton, 
hito adapisoJocal conditions by taking noo-statonary 
mean and variance ol tre image, The non-statignary mean 
escribes Ino gross structural the imag, Wii varianco 
provides fst order edgs and texture information. In ts 
context and based on Wiener fitring (Lim, 1990) Iramo- 
work, a formula is derived that has tho same: form ana! 
swucture ol the original Loe's muttiplicaivo but with var 
ab variance nos a (), his fitor was applied to both 
Images violin he best visual results, 

lis important to point out tha ho valuo lor ¿comes 
rompo taken direct from the SAR images. tmay not bo. 
he same lor serial or satolta images; each class requires 
hs cin estimation ol o. Nowse speckle sn general oltho. 
order ol magntude ol ne detal, so increasing a ncro- 
ments tha degroe o! fltring but emears fino detal 9 ol, 
Semilan, the densty function associated lo Ino spacklo 
"ola is direc related to th number of looks. 

Inia fnalstop. speckle reduction can bo carried 
ut as well by means of a geomerri fiter based upon te 
Hitor Miss (Prat, 1991) morpnological operation. Several 
sio assumetions are made lor this ter, First ol al. ho 
úspecile nolse pattern 1 hought o be Iormed o apikes ol 
no more than 3 pixels wide. This isin agreamant o the fact 
trat the autocorrelaton function or specklo tados out very 
rapidiy lor second or farther away neighibors. Obvious, 
getals ol 3or less pixels wil be wiped out. Second, he 
image fed is supposed o be staicnary. And thirdly, the 
noise and the image are uncorrelatod since the spocklo a 
Independenty generated lor each pixel, The Hit or Miss 
transtormatíon is defined (Prat. 1981) as ollas: A binary 
valued pattem mask is constructed over a given window. 
“As the window sides over the image, Ihelogial stato (Dor 
1)oftba image pixel in corespondenca to 1h center of Ino 
éndow 's changed whenever the pattern mask matches 
he pixel alves in tho image. For he morphological ter 
the llowng eight masks are considered 
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“The geomerrc fiter algorihm works as follows: four 
¡directions are defined, horizontal, veical and two main 
'Siagonais, along each one a gray value prole is consid- 
eredas a binary image: the ares under the profile, defined 
as the umbral, is taken as ones, the restas zeros, To th 
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binary image, the Hi or Miss transformation is applied 
'employing tha masks above definea, thenthe complement 
of Ihe results obtained and the masks are applied agan. 
"Note that oniy a certain pattorm ol pieels is axplcaty given 
lor each mask, Whenever this pañtern is encountered in the 
umbral, Ine HitorMistranstormationis applied regarciess. 
al he value of the rest of the porels. Such procedure la 
úxeculed to every prolie foreach girectionandthisistaken 
¡as ono fteration of the algoritnm, Several terations are 
required to elminate the speckle, and a quanttaie crte- 
ion is provided to determine the-.end of the algoritwn. 
vonthough tha procedure is convergent. On tha other 
hand, this partcular Hit or Miss transtormation may be: 
Viewed as an specíf case ol erosion and ailon (Prat. 
1991) tnat wipes out the spikes o! the umbral and hence: 
réduces Ine speckls. This fiter provides Ine best visual 
results (Figs. 3 and 4), so ltwas selected as the previous. 
step to texture segmentation. 

Lotus asume now that acortan homogeneous area in 
he imago contains on|y one object, For area a certain 
mean yu and standard deviaton a aro associated. As the 
ú'spocklelareduced ineachiteration, a lower value ot amay 
be expected, so the ratio a= aja sa rs degree indicator 
le smocihness ol Ine object area: is zero fora Nat area 
andtakos a certain maximum valve tor an uníitered caso, 
Foreachiteration tne spechle indexa s evaluated, when ts 
chango value lor succosvo lterations ía loss than somo 
predelinad Ihreshold, 4 «$ e ne algoritm is terminara. 
he speckloinciex may be defined as well as an average ol 
rios fa tor non-overlapping windows covering the im- 
ge. ILappears hata estimated from homogeneous area 
provides Iho best measuroment lor ino degree ol specióo 
reduction, eventhouqh It requires a manual Reratin at the 
begíning ol the algoritm. Let us point out final nat the 
¡onceptotumbral sony lo lustration purposes, in practica 
ho algorilm works lrectiy wit image pica! values along 
tha our mentioned directions. 


4 - SEGMENTATION BY REGION GROWING 


Speckle reduction n the present context is ín fact a 
proprocessing step towards. segmantation ol tha image 
Into comprehensive non-overlanping. regions. Specie 
noise cificultne classfication ol mago pels in meaning- 
Ibi categories, hence ts modeling and fitering. Its won 
tomentionhalas edge preserving itecng isapelecto te 
image the noisa is reduced and ha histogram translokr, in 
¿eneral, rom unimodal to multimodal preparing Ine image. 
lo better segmentaton. Segmentaton ol an image is ne: 
grouping of pixels ínto connected and non-overlapping. 
regions thal may possess a meaninglulcorespondence to 
ú>bjectsin Ine scene. The grouping 's carried out based on 
spectral and spatial propenies of pixels. The output e 
gions may not have a clrect or evident relation to scene. 
bjects represented in the image. Segmentaton sn any 
case a necesary step towards image understanding and 
pattern recognition. 

Region growing (Bes! 1988) s one ol the methods lor 
image segmentation, Many radar images are monoctro- 
'matic $0 they are suited for region growing techniques; a 
lormaldefiitonnow allows: Given animage tana alogical 
“uniform predicate Ele) ina parition or segmentaton S 
of the image in terms of a set of regions Al Let N, be the 
number ol regions in the segmentaton, Then the segmen- 
tation conditions or the set Sis established as 
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1 qt ote Lv w 
lat 

2-ANR=0. ve * a 

3- Asa 4-comected set ol ports. 0) 

4- AR) 5) 

S- ElRUAj=false, Wim). 6 


“he result ot he segmentation s he set S given by 


N 
Ss UA (6) 
l=1 


The key aspect in this region growing schemo ls 
ú"evicenty tha specific delito al the unormitypredicato 
Elelandsomaconstrains mposedon he process Tha ir 
restrictions 1o require hat Als composed ol a minimum 
rumber ol pixels order to avokd tha Invi!satistactión of 
precicate (4), Other restrictions may include a certain 
shape or maxmum izo lor A The growing ol region may 
stop as well when some ego levels encountered. Evi- 
¡enty ho more restrictions, the higher ho computational 
¡cost Forconstalnedlow-lvol (Bes), 1988) digital analysis 
systems, segmentation can be achieved quick and accu- 
roto, but for less constrained envirorments, high level 
semantic knowledge about a scene ts usualy roquir in 
order 1o oblain a precise segmentaton wi meaninglul 
image regions 
Three crteriawere adopto or undormiy predicato Ele 
1 Tone Ditterence - Tho ciforonce valuo, betwoen a 
now contered pixolasna tha averaga olitsmolgnbors 
isevalunts, the result does notexcosda predofined 
"broShold, the conter pines agragated to the region 
2- Correlation - The cfferenco value botweon a window 
entered pixel and the normalized correlation between 
ire onginal ana tho actual Windows a evaluato, tha 
result does nat exceed a predefned throshold, the 
¡enter pixel ls agregated to tho region 
3- Entropy - Tho diferenco entropy value botweon tho 
original and the actual windows s evaluatod, ll Ino 
result does not exceed a predeined Ihreshola, the 





Center piel 's agregated to the region 
[na coectancn api oa ogin A depa 
on the closeness of value o given untormity ertrión 


Inany case, the growing commences by seeding. pixel ln 
a homogeneous area and stops untl the unitomity proc 
ale isroJonger satis, The algoritm works tor variable 
windows and predefined thresholds. is wr to mention 
that SANDI (Lira and Landgrave, 1992) was used as a 
based environment to instal and test the methods deve 
ope in this work 
In the present stage ol development he region grow. 
ing procedure begins when a very typical pixel ola glven 
ype object or region is seected in the image. The position 
ol seed pirelsarenoterucialas long as the Ihreshola factor 
is adjustod according o the mean and standard devíaion 
of the seeded window neighbomood, As the window size 
is increased coarser textures may be taken into account 
butt the expense of reduced cerainty as 1o the position 
and definition ol re classes. Increased ceralnty ol class 
memberstup is at the expense ol spatial accuracy. Class 
fication uncertainty s tilo be worked out 
“Thvee aspects ere crítica inregion growing; tho seed- 
ingota pol, hesizeol the window and the threshola. The. 
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firtonels done intaractvely by selecting a representatwe 
plsel ol he class texture to be segmented. Centered on 
'sach seeded phael a 1 x 1 size window ls considered lor 
which e! /4) is a fist approximation 0! dispersion pixel 
values. Only odd windows are employed so | = 2a + 1 
'uhero a is natural number. The parameters a ana y) aro 
he standard deviatcn and he mean of he window piel 
Population. The location of a sesded pines not critical as 
longas the Ihreshold is acjusted accorging to na Jollawing 


e, Kal, forcloss m 


vwero K, ls.a consiant that normalzes the tvreshold in a 
range belwsen Ong 1, Inother words, fora given window, 
¡La region growing is launched from Iwo dferent porels 
belonging!othe sama tone class (Tone Dillerence).a basic 
equivalent segmentation is obtained provided the trrest- 
¿ld is ajusted1o local winciow conditions according 1o (6). 
ha range for Ine hreahold is between O and 1, zero. 
value leads to Ino tval case ol an empiy class, the one 
value yields (ho ente imago; a best estimation for ria 
vroshokd Wi be discussed later. Asa ft instance a best 
erivallon for window size is presented for Irae unitormiy 
precicatos employed in Invs work. Lat us constuet a 
Iroquency distribution pixel values (Gong and Howarn, 
1992) for ch iniial window class ol seeded: pixels. Let 
A de a clistribution for class | and window size 1x1 
where v is a pixel volve and the function Mv) sell ts 
Iroquenoy ol occurrenca, A distance between a pair (1,91 
classes may now bo delined 


$ 
dp Sy hh ws 
a mr 

wo 


where N,ístho highest possible pixel value. in Ihe prosent 
casoN, =255.The average distance lor alpalrs ol classes 
istnon 


Sa 
de boe h (9) 
el dew 
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hora, ista number f classes. Fortha Tor Dilorance 
Uniformiy precicato, he best window size is given or d 
maximum. 

The Correlalion an Entropy unilomity predicates. 
special suled or texture analysis, deserve a somenat 
lterent approach. («JJ and (min) are te cooreínates ol 
seoded písels or classes land rospectvoly, ho corrola- 
lion between window classes 2, and 2, ls 
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ho averago correlation for all palrs ol classes is then, 
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he best window size is Ine one that generates ne 
mínim €. 
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The Entropy fora given window class is 


ey Hee, o) 
SS 


where C, are the elements ol Ino co-occurence matrix 
averaged over the elght possible directions; a detailed 
account! trisformula ls given inthenext section. Tololow 
he same ratonale as above its required to establish an 
entropy value for a: pair of window classes. Tha 
¡co-sccurence matrix s generatedon the normalized num- 
bero! pixel pairs that appear in a given crection, Far two 
classes let construct a co-occurence matrix based on 
Pals thatoccut pixel by pare on ho wo1windows, For two 
dentical classes, tnis co-occurrence matrix would: havo 
only onenan-zaro element equaltoone; the jointentropy! 

would be zero; lor diferent classes a contain valuo wow 
be obtamed reaching a maximum when th classes aro 
total diferent. Theratora,tho avorago value 


e El (13) 
AN 


NS] 
+ maximum for ho best window sizo, 





“An initial ostime ol the throshold: Comos from each 
uniormty predicale and from local conditions ol seedad 
plaels windows, terefora or a certain class window with 
Center coordinatos (k,), tve correlation Ireshold is 


A ME 0) 
9% 


ad lor Entropy 1 
AI (15) 
mos 


The throstoká lor Tone Dilrence ls already given in 
'quation (6). The constantsk, nah, re or norrmalzation 
purposes as discussed above. Tho best Ireshold is de- 
ived alter examirig the Jetries Malu, matrix tor al 
stances among classes, The generationo! atexturamap. 
equis a giobal categorization ol texures for the whole 
Image. 1 was done In two stops: an ontropy evaluation, 
anda clusterig analysisovera hybrid mutispectral mago 
composedtoliwo band: the crigial image and the output 
fine entropy analysis. 


5. TEXTURE SEGMENTATION 
OF ENTIRE IMAGE 


Astandalone texture segmentalion was developed by 
means ol he entropy crterion (Haralick, 1979) based on 
the co-occurence matrix (Gonzalez and Wintz, 1887) 
concept, ín combination with lustering analysis. The 
co-occurrence matrixis constructed as falows; Let Bbe a 
vector position operator letk be Ihe number o alowable 
plselvalues inthe image, and lt 4be-ak xk matrixwose: 
(6.2) element s the number of times hata plrel having a 
Value A, occur a position given by Brelate to a pixel 
having 3 value 2, The operator Bis defined here or the 
sgrt directions Gwen by the Freeman Code. Since the 
correlation between pixels fai of rapidly as hey become 
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apart, 1ha longih 'oT'£ noed not 10 be greater than: tha 
distance lor closests noighboes. Thus, lo given pixel, he 
B-comected neighbors usualy contribute he most to the 
mato Note hat tho dimension o ys doterminedd by the 
iraxmum pieel value regardiess ol image size, In the 
prosor implemantaión a lcing wincow ol variable siza 
(9x3,5%5,...11x 11) may be employed to calculato ne 
hi, matrix associated 10 cach pixol neIghiborood. The 
element mol stores tho numborof times nata pair ol 
Values (2, A.) oCCur In relate position given by once la 
ls calculated, each element m,, ls divided by he total 
number of palta, tv resul 1 a noWmatrx Came Iho co- 
occurrénce matrix, Ihat characterizes in 00d 
opproxmation lhe jolt probably that a pair o! pixel 
salisiying Ohavevaluos (A). ln olher words Ino marx C 
1s an approximation to 1h olnt density function, There ls 
axial Symmetry in Iho eight directions chosen, so Cs a 
symnetio matrix, This manes 1 facto bidimensional 
histogram when tho window equals tho image szo. 

To characteriza the content ol the co-occurronos m 
tu, ho ONIsopy crteion was ulizad 








lo E 
Envopy= - E, Y G.LOg5, 
pa 


were 6, aro ho alements ol matrix C. Tris expression le 
based on Iharmociynamics laws and provides a moasure- 
iman of ofder and probablity ol a closed systern ín ho 
Following way 


High Entropy» Las Order Higher Probabty ol Oecurronco, 
Lowor Entropy - More Ordar » Lowwr Probabily ol Occurrenco, 


Te Jr denaly function la alwayt loss inan or 
qualto'ono, Wen piro voluas to alcing window aro 
“leual, 9ne of NG, = , and toros ase zoo, 
Entropy 0 loweni poseo valuo no pix val 
all iforen, non, < 1 and Entropy > O, wi a cenaln 
rraximum value, Therelora for a, smooth textura (high, 
¿ror) oa no imag, low pixel aluo zona (dono) 
correspond in tno texure imago, A coarser extur tr 
luces into cloar tones, this Implias What a toxtur 
sogmoniation map ls posible by clsterng no tonos (pix! 
values)ofIhetoxtue image, This was: doneby means ol re 
tollwing hypottosis: ho original imago contains ho spec- 
tral (one) information o re scene, and tha texure imago: 
caries tho spatial information Iereforo a combination ol 
ihesa Iwo Imagos Into a 1wO band muliapeciral Image 
should provide the texture analysis ol Ihe scene, Tho 
analyeis was carried out win on Isodata clustring algo- 
ttm 




















6. RESULTS 


TWO SAR Images ol 512x542 pixels, vit resolution 
dl 25m waro usod 10 lost the melhods developed in his 
work, Ones a subreglon o! a SEASAT image known as ho 
'Obayos anticlinelocalegin norinern México, the olher one 
is an aerial image of norther Holland. (was contimad 
experimental tnat the speckle density function lor he 
'SEASAT image ls gaussian, and Raylelgh ta he aerialone, 
To both images a serios ol fllers were applied (section 3) 
wihthe above discussed results. Eventhougn Lee! Mli- 
picalive blures somewhal the edges of 1ho image, lt 
performs as good as Geometre fiter as long as region 
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¿growing is cencermed. For entiapy texure classiication 
Plxposes s entero uso he goomatric ler n acvance: 

Twolypeso! region growing procedures were applied 
t0bolhimages,Forihe SEASAT image, Iwo textura classes 
vere cleary distinguised using Ihanormalized correlation 
orion and a 5 x 5 window, such classes presumably 
correspond to two diferent Inologíe units, ore wr 1 
this drectón tsreguirea. For the aeril image (Flg.5)a oty 
anda cultvated field oro clgary separated employig o. 
same crtrion and same window, Some tests Wih 1h 
+ntropy crion were also performed anda vinciowo! 5 
$5 sooma to provide a moro delalled. sogmemiaton for 
lexturo oventhough he Computational cos! ls higher Al 
tos Wero applied after ono Heration l he, geometria Mer 

A separate classification by means ol co-occurrenco 
sti ana entropy csterión was worked ut or the wal 
SEASAT images A subsequent clustering operation wi 4 
elustersindicates a cioar categorizaion Fig 6)oftextures. 
ho bands amployed in this clusterino wero he tere 
maga by gaomatric filter wi on Neratlr (a= 10%) ana 
ho one resulting rom entropy measurement. A. compar 
¡son'was done wit clustarina (Fig, 7) ullling the origina! 
¡imag wit no fítorng and lhe one resulting from entropy. 
Monsurement In both casas ho covárlance matrbcanid the: 
¿Jelfios.Matusita istanca (Tables 1and2)were evaluated 
A Visual inspectión (Figs, 6 and 7) and a. comparison ol 
Ihoso matricos show a claar improvement ot lasificaion 
hen the fMiored image is ncluced In the clustering analy- 
sl 


7. CONCLUSIONS 

















A meinod to detact-spooifc toxtura:objecta. or to 
¡general an overalltexture segmentation ón SAR images ls 
prosonted. A rationale for specke reduction andits impact 
on toxturo analysis is givon as well. Rosulá aro prasintad. 
armploying and seal and a satóllto SAR images. buho. 
method Is applicable to any rango Image goneratod by 
means ol coherant fadiation 

Incomparison o pravidus works (Corners, Thvedi and 
Harlow, 1984; Pedi and Frarkin, 1891: Ho and Wang. 
1990), ino ntropy calculations here aro based on the co- 
Gccutrence matrix averaged over tha, olght diraclions. 
defined by tha Freeman Code, This reduce Io lexture: 
dependence on possible variation ol luminance cond:- 
"on or Other monotonical sti ol pixel valvas over tha. 
image. 

Detectloncltextureobjectalizing theragin growing 
procedure may bó seen us a previous stop lo evaluate 
morphological and topological characteristics ol selected 
'aspoctsolthe scene; These data mayin ur be usadas an. 
Input o a dalabase lor management ol resources related! 
to tho scene. 

Reduction of speckle has a:clear consaquence on 
texture segmentatión: here 1s a visual Improvement alter 
(tering. resulting, ia /amoot map. (Fig. 6) ol toxtures 
a ragmentanon ol (lg. 7) ol 

Inacddlloo to is, Ihe covart- 
“arce mat and. ho Jeffries-Malusita distanco show an 
improvement n compactness and separation of classes, 

Ageneralizationolentropy evaluationtomuliespectral 
Images could provide a contextual:spectral class how. 
ver the calculations would be quite involved yielding a 
high Computational cost 
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TABLE 1.- COVARIANCE MATRIX AND JEFFRIES - MATUSITA DISTANCE BEFORE SPECKLE REDUCTION, 





ase Texte Covaranoe Mata tros Matsta Diane 
EEES e 
ERE 
oo ora om oo 
Four (o) ZE 288 oO O 10 o 
Intermediate (3) — 258.04 27.33 A a Bio der Bit 
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TABLE 2: COVARIANCE MATRIX AND JEFFRIES - MATUSITA DISTANCE | 





Smooth (1) 1789 
4822 14075 
Regular (2) 11590 assgr 1 0.000 101520002000 
24427 — 188023 2 1015 0000 19182000 
Intermeciate (3) — 49853 25225 3 2000 1918 0000 1975 
25225 47864 4 2000 2000 195 000 
Rough (4) 1752.33 — 125493 
1125433 154213 
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Resumen 


Las Imágenes de radar de apertura sintética son en general monoespoctralos 
proporcionando un rico contenido de texturas de la escena, Debido a la radiación 
olectromagnética coherente empleada en la generación de imágenes SAR, un ruldo 
;pecklo") se sobrelmpone a ellas oscureciendo la Información textural. 








En la prosente investigación, se presenta un método para detectar (segmentar) objetos 
'quema de crecimiento de regiones con predicados de 


toxturalos empleando un 
uniformidad especificamente 
toxturas de toda la Imagen, se propone un método hibrido: se aplica un criterlo de entropía, 
basado en la matriz de co-ocurrencia, creando una banda sintótica que os incorporada la 
original, on un conjunto multiospectral de dos bandas, Se roaliza un análisis do cúmulos con 
lun Algoritmo de Isodata sobre este conjunto generando un mapa de texturas. Una Imagen 

iálto son utilizadas para ejemplificar la metodología 














So proporciona una discusión sobre la reduccion del "speckle" y se proporcionan 
detalles sobre su Impacto ena segmentación. Ellrajeal"speckle” e pllcnalas imágenes 





tanto antos como después del análisis de toxtun 








1. INTRODUCCION 


as imágenes de radar de aportura sintética. 
(SAR) prosentan un estos dias un gran 
Potencial de aplicaciones en muchas áreas de 
la Invostigación (Kolawoio, 1991; Kwok eLal. 
1991; Marom y Shemox, 1991; Loberl, 1990, Wang y Ha. 
1990). El contanido de información de las imagenas SAR 
nplementaro al proporcionado por los sensores. 
OplicoscomoLANDSAT y SPOT, Siendoln sistema activo, 
el SAR es capaz de luminar el área del blanco bojo 
dilorontes ángulos de vista, varios esquemas de 
polarización y longitudes de onda selectas, La amplitud de: 
la radiación rotiocispersada depende básicamente de la 
rugosidad y pendiente delblanco, delangulo:ncidante, de 
la Songilud de onda y de la constante voluméirica 
¿olécirica la que, asu vez, es función dé las propiedades. 
físicas del blanco, Las imágenes SAR muestran un gran 
detalle acerca de la superticie del blanco, oscurecido 
Unicamente pot el efecto conocido coma "speckla 

Las Imágenes SAR ganeran un ren de ondas de 
pulsos electromagneticos coherentes y en fases que 
interaccionan con la superficia bajo estudio, Las 
caractorisicas geoméliicas del tren de ondas (Elachi, 
1988) son tales que, unaceldila untariao de resolución es 
lurinada por la radiación y la Intensidad dal pulso 

















Trabajo apoyado por al Proyecto INIOOZO2 de 


retrogispersado determina ul valor del plxel asociado ala 
¿olclla, No obelanta, enla celdila oxiste un cierto número 
de objetos con formas y orientaciones aleatorias que 
“aCIuaN COMO dISporsoros quando sx lamano as del mismo. 
ren de magnitud de la longitud de onda del tren de 
ondas incidente, Estos objalos retradiapersan ol ren de 
ondas Incidente en ciferentes areccionas y fases. Como. 
Unrosultado de osto, ocurre un fenómeno ue inortorencia 
que genera interacción destructiva y constructiva de la 
radiación, El fonomeno de interacción us conocido como 
lacio coherento de "speckle”, smilar en apariencia en la 
imagón, al grano de una emulsión fotográfica, Como una 
delinición, al 'spackle" (Dawaolo, el al, 1989) es ol 
rosultada dela Intorferancia destructiva o constructiva de: 
las ondoletas cohorentos destasadas ontre allas, 
rovodisparsadas desde una superficio cuya rugosidad ws. 
comparable a la longitud de onda de la radiación 
empleada, Como tal, el "specke" es dependiente de la 
textura local de la escena bajo ostudio. En esto articulo se 
presentan una serie de resullados para segmentar 
texturas en imágenes SAR, empleando un esquema de: 
crecimiento dle regjones para la delección de objetos 
texturalesespacilicos y una medida de entropía basadaon 
la matrz de co-ocunencia, para generar un mapa de 
texturas itoctamente relacionado a la superticia de la 
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2. REDUCCION DEL RUIDO "SPECKLE" 
EN IMAGENES SAR 


Una manera de reducir fisicamente el "speckle" y el 
ruido aleato es porel promedio incoherente de N.vistas 
IN =4 para el SEASAT) de la misma escena; esto reduce 
o obstante la resolución azimutal de la imagen. El 
promedio incoherente es el cálculo dela media arimérica 
de Nistas; este tipo de promsdio borra la Imagen y 
r6duce ol detalle y requiero de varias observaciones dela 
escena que no siempre están disponibles 
Laotra alternativa pará la supresión dol'speckle*es 
por ltraje dígtal de la imagen: esto demanda algunas 
suposiciones en cuanto al comportamiento 
;peckle". Para imágenes de Intonsicad de una vista la 
unción de densidad del campo del "specklo" está. 
descrita por una función exponencial negaliva. Cuando. 
so consideran imágonos de Intensidad de Nvistas 
(Leberl, 1990) la unción e densidad es más bien una 
Chiroualdrada con 2N grados de Ibon, a medida que 
N se Incrementa, esta función se aproxima a una 
distribución normal. Teóricamente (Elachi, 1988), el 
ruido 'spocke* para imágones de una vista obedece a 
Una exponencial negativa y de aquí que sea 
Muliplcatvo en el sentido que su desviación estándar 
8 proporcionala suecia (Loo, 1983) Noobstanto, a 
sido argumentado (Frost y Yurovsky, 1985; Tur el al, 
1982) que el modelo múltiplicatvo no slempro os 
correcto. Por otro lado, el "speckla" no slompra está 
considerado como totalmente desarrallado (Nezry et 
Al, 1991), ya que para áreas urbans en Imágenes de 
ala resolución, solamente algunos dispersores están 
Inoluldos en la cell unitaria 
En particular, las Imágenos SEASAT y ERS=1 son 
promedios de N.istas yrosullados experimentales (Les. 
1081) muestran que la función da densidad os gaussiana, 
Resultados experimentales (Dewaele el al. 1989) 
Muestran tarmbión que on imágenes SAR, los pixels tienen 
Una paqueña correlación:con sus mas cercanos vecinos y 
prácticamente se desvanece para los segundos yerceros 
nos; como resultado le esta, ruido "speckle" puede 
"considerado como blanco. Aun cuando el "spackle" 
está tocricamente muy corrlacionado, en la práctica la 
suposicion de ruido casí blanco para imagenes del 
SEASAT y ERS-1 puedo ser el resultado del promedio do 
'Navstas Silacorolación lapa aser despraciablo, sebo. 
giseñar un algoritmo de corrección con media y varianza. 
locales y dado que el "speckle" está directamente 
ralacionado con la textura lacal ol algortmo deba ser 
AaGapIaDlo 
Para un 'spechlo"no correlacionado con exponencial 
negativo como función de densidad un modelo de 
corfección Iuo desarrollado (Kuan, 1987) que conduce al 
llamado fitro máximo a-posteriori (MAP), Otro modelo, 
basado en datos experimentales (LI, 1988), sigue una 
distribución gausslana no-blanca razonablemente tien al 
histograma del "speckle". Basado en resultados de 
correlación entre pixels vecinos, se adoptó un modelo de 
*pecie" dependientemente correlacionado con a señal. 
El tamaño de la correlación cel "apeckle", es en este. 
modelo, aproximadamente del tamano de la función 
incoherento de dispersión de punto del sístema. 
Los modelos del "epecke", hacen de una manera u 
¿tra, varias simpilicaciones (Leber!, 1990: Dewaele ot al. 
1989) ue en la práctica no funcionan, 
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La onda: electromagnética. Incidente. está 
perlectamente polarizada. 
2» La onda electromagnética incidente es perfectamente: 

monocramática. + 
3. Él número de dispersores en la: cell 

necesariamente grande. 

4» Elnúmero.de dispersores en una celda unitaria: 
constante. 

5. Las fasas asociadas a los dispersores tienen una. 
unción de densidad unilorme, 

Esto significa que los fos basados an modelos de 
"speckle' no se comporian tan bien a medida que se 
desvían de estos hechos. Además de esto, los fitros 
Venden a borrarlos bordes de la imagen ya que siempre 
introducen algun tipo de fraje pasa bajas. El objetivo es. 
reducir el. ruido, preservando en mayor grado: los 
scalones, los hordas y las Inaas de la imagen, Ellipo de 
tro para esto, depende de lalunción de densidad o dela 
ú'struclura geométrica del "speckle"; más delallos de ésto 
o lan en la siguiente sección, 


3. FILTRAJE DIGITAL DEL "SPECKLE" 


untar os 














Comenzando: con un enfoque simpla, vanos fltros 
acioptablos fusron aplicados a dos Imágenes (Figs. 1 y 2) 
SAR solacias; aorotramsportada y satoltal Suponiendo. 
que en regiones homogéneas la población de plsels 
Presenta una distribución de valores aleatoria unllorme y 
¡que en ároas do bordes la población de pixels pertenece: 
en su mayra al, se aplico un promedio restringido y. 
ponderado sobra una vontana, esto ue ol caso para los 
fltsos de k-próximos vecinos y el sigma, Habiendo 
Aplicado estos fttos para ventanas do 5x5, pormanaco un 
lo dl ipo moteado, por lo que una cascada con fito 
de mediana parece apropiado. La combinación de 
cascada Sigma + Mediano, ara ventanas de 5x5y3x3, 
respectivamente, produce resultados razonablemente 
buenos, aun cuando los bordos son un poco borrosos, 

En un segundo paso y suponiendo que el "speckle" es 
na lunción do densidad gaussana no correlacionada y un 
ruido detipomuliplicalvo, estoconduce ala aplicación del 
Állro e Los mullplicatvo (Leo, 1981) con una ventana de 
5x5 y una varianza del uldo de o =0,28, cuyo valor está 
tomado a partr do datos exporimontalos. as dcir de la 
Imagen misma. Resultados de esto fro muestran una 
pequeña mejoría visual con respecto a la cascada 
anteriormente méncionada, Con un modelo de "speckle" 
cependienteménte corelacionado con lasenal so produce 
Una generalización del fito de Leo muliplicalwo. En esta 
generalización tro e adapta a las condiblones locales 
tomando una media y una vavanza no-outacionarias de la 
imagén. La media no-estaciónaria describe la estructura 
gruesa dela Imagen, mientras que la varianza proporciona 
información de primer orden de los bordes yla textura. En 
ste contexto y basado en el marco del itraje cie Wiener 
(Lim 1990), se obtiene una fórmula que tene la misma, 
dstructura y lorma que lacrigial del Les multpiicativo pero 
¿que tiene una varianza variable para elruido. 11), esto 
flo se aplicó a ambas imágenes proporcionando los 
mejores resultados visuales, 

Es importan puntualizar que el valor de ¿proviene 
degráficas obtenidas directamente delas imágenes SAR, 
puede no ser el mismo para imágenes asrotransporadas 
o satoltales, cada clase requiere su propia estimación de 
a; El rudo "speckle es, en general, del mismo orden de 
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magi que los detalles nos, de tal manera que al 
Incrementar. 9 se incrementa el grado de fitrae pero 
también se borran los detales fos, Smiamente. la 
unción de densidad asociada al campo del"specide" está 
directamente relacionada con el número de vistas 

Ent tercero y úlimo paso laregucción del"specide” 
puede lavarse a cabo también por medio del fito 
geométrico basado en la operación morológica (Pan. 
1991) do Intersecta o No-intergoca. Varias suposiciones. 
básicas se hacen para este lo. Primero, el patrón del 
1ido"speckle" se piensa como formado por picos de no 
más de 3 plot de ancho. Esto está de acuerdo con el 
hocho de que la unción de autocorelación del "speckte" 
se desvanece rápidamente para segundos o más alejados 
vecinos. Obviamente, los detalles de 3 0 menos phels 
serán eliminados, Segundo, el campo de la imagen se 
supone estacionario, Torcor, el ruido y la imagen estan 
descorrelacionados ya que el "speckle" se genera 
independientemante para cada pixel. La transformación 
Intorsecta o No-ntersecta se celine (Pra, 1991) como 
gue: So construyo una máscara pato binaria sobre una 
vontana dada, Ámecida que la ventana sa desliza sobre 
laimagon, olostacológIco(06 1) del piel de a magenen 
correspondencia con el centro de la ventana se cambia 
todavez quel patrón máscaracomcide conlos alores de 
la imagen, Para el lvo merológico se consideraron las 
siguientes ocho máscaras. 
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El algara dello geométrico trabaja como sigue: 
se definen cuatro acciones, horizontal. vertical y las 
lagonalos principales, alo argo de cada una de ells el 
porfi de valores de grises se considera como una Imagen 
binaria; ol droa bajo ol pri, dotido como umbral, se 
loma como Unos, el resto como ceros. A. esta Imagen 
binaria, so lo aplica la translormación Inersacta o No- 
intorsacta empleando las máscaras arriba. dofinidas, 
entonces se obtiane el Complemento del resultado y las 
máscaras so aplican nuevamente. Nótese que solamente 
un cet patrón de pixols os axpliotamente dado para 
cada máscara, toda vez que esto patrón se encuentra en 
+l umbral, la translormacion Intrsecta o No.intersecta es 
aplicada independientemente. del valor del resta de los 
lts. Tel procedimiento es ejecutado a todo perti para 
cada cirección y sto es tomado comuna interacción del 
aigorimo, Varias raciones son requeridas para eliminar 
el'specido'yuncritriacuantitaivo.es proporcionado para. 
determinar el. fin del algoritmo, aún cuando 
procedimiento es convergente. Por oo lado, la 
lranstormación particular Intersecta o No-Itersacta puede 
ser vita como un caso especifico de erosión y distación 
(Pratt, 1981) que olimina los picos del umbral y eciuco por 
lo tanto el "specid" Este tro proporciona los mejores. 
rosulados visuales (Figs. 3 y 4), de tal foma que lus 
seleccionado como un paso previo a la segmentación 
textura, 








Supongamos ahora que una cierta área homogénea 
en la imagen contiene Unicamente un cbjeto. Para esta 
Area una cierta medía y y una desviación estándar erson 





manera que el cociente a=0/p es un indicador de primer 
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¡grado dela suavidad del área del objeto; es coro parauna 
área plana y toma un cierto valor máximo para ln caso no 
litrado. Para cada iteración el indice. del “speckle" a es 
evaluado, cuando su valor de cambio para lteraciones 
sucesivas es menor que un umbral predelinido a as e el 
algoritmo es terminado. El indico de "speckde” puede sor 
¡definida también como un promedio de cocientes 0/4 
para ventanas no traslapadas que cubran la imagen. Todo 
parece indicar que a estimada a parir de una 
homogénea proporciona la mejor medida del grado de 
reducción del"specide", aun cuando esto requiere de una 
iteración manual-al comienzo del algoritmo, Pun- 
ualzamos, Ínalmente, que el concepto de umbral es 
únicamente para proposttos de lustración, enla práctica 
el algorimo trabaja directamente con los valores de los 
"purels de la imagen a lo largo de las cuatro direcciones 
mencionadas. 





4. SEGMENTACION POR CRECIMIENTO 
DE REGIONES 


La reducción del "speciie" en el presente contexto es. 
¡de hecho una etapa de proproceso hacia la segmentación 
¿e la imagen en regiones significativas no traslapadas. El 
Udo 'specide" aficula la clasificación de los pixels de la 
megenen categorías con significado, de aqui su modelaje 
y rado. Vale la pena mencionar que al aplicar un tajo 
ue preserve los bordes en a imagen, el ruido se reduce: 
y ol histograma se transforma, en general, de unimodal a 
Multimodal preparando la imagen para una mejor 
segmentación. La segmentación de una imagen es el 
agrupamiento de los pixels an regiones conectadas. 
úo=raslapadas que pueden poseer una correspondencia 
significativa conos objetos dela.escena, El agrupamiento 
so lleva a cabo. basado en propiedades espectrales y 
espaciales delos pixels. Las regiones resullantes puecion 
únotoner una relación drecta o evidente con los objetos de 
la escena representados en la imagen. La segmentación 
ds, en cualquier caso, un paso necesario hacia el análisis 
e la imagen y el reconocimiento de patronos. 

El crecimiento de regiones (Besl, 1988) es uno delos. 
métodos para la segmentación de la imagen. Muchas. 
imágenes de radar son monacromáticas por lo que son 
adecuadas pera las técnicas de crecimiento de regiones; 
una definición tormal sigue añora: Dada una imagen! y un 
predicado lógico de unilormidad Elel. encontrar una 
úparticióno segmentación Sdela imagen en términos de un 
conjunto de regiones A, Sea N, el número de regiones en 
la segmentación. Entonces las condiciones de la 
segmentación para el conjunto S se establecen como 





N 
1 ÚR=1, donde REL. YI m7 
fal 


2- ANA=0,Viaj a 
3- Res un conjunto de pleels 4-conectados. — (3) 


4= ER] = Verdadero, YI 0) 





ERUA)= Falso, vi) 6 


SELPER 
El resultado de la segmentación es el conjunto Sdado por 


N 
s=UR 





El aspecto principal en ste esquema de crecimiento 
¿e regiones es evidentemente la definición especica del 
predicado de unilormidad Ele) y algunas restricciones 
impuestas al proceso. La primera restricción es requerr 
que Resté compuesta de un número mínimo ce pixels con 
objeto de aviar la satisfacción vil del predicado (4) 
Otras restricciones pueden incluir una cierta forma o un 
tamaño máximo para A, El crecimiento de una región 
puede parar también cuando algún nivel de bordo sea. 
encontrado. Evidentemente, entre más restricciones, 
"mayor el costo computacional: Para sistemas de anaiais 
digital de bajo nivel de restricción (Ses). 1988), la 
"segmentación puede alcanzarse rápida y precisamento, 

medios: manos. restringidos, se requiere 
un conocimiento semántico de alto nivel 

Acerca de la escena con el objeto de obtener una 

segmentación precisa con regiones significarvas. 

Tres crios se adoptaron paro el precicado de 
nllormicad Ele] 

1 Diferencia de Tono El valor dierercia este el porel 
centrado en una ventana yl promegiode sus vecinos, 
es evaluado, si ol resultado no excede un umbral 
predefido, el pbrel central es agregado ala región. 

2- Correlación=Fivalorcioranciaontroal pol contrado 
+n una ventana yla conelación normalizada entre la 
úvontanoriinal yla actual es avalado, sí elresultado 
no excede un umbral predoliido, el piro! contra! 
agregado a la región 

3- Entropía - El valor diferencia de entropía entra la 
ventana original y la actual es evaluado, elresulado 
o excede un Umbral predefinido. el pel central es 
úagrogado ala región 
En brave, la aceptación de un pixel una región A. 

doponde de la corcanía de su valor a un crtero de 

uniformidad dado, En cualquier Caso, el crecmiento 

Comienza al sembrar un piel en una rea homogénea y 

pora cuando el predicado de unlormidad no so satslaca 

más. El algorimo trabaja para ventanas variables y 

Umbralosprecéfinidos, Vale la pena mencionar que SANO! 

(Lia y Landgravo, 1992) fue utlizado como un medio base. 

para Instalar y probar los métodos desarrolados en el 

prosento trabajo. 

Enel presente estado de desarrol, el procedimiento 
dh crecimiento de regiones comienza cuando un pixel 
pico de un objeto oregión tipos sembrada ena magen 
La posición del piel sembrado no es crucial mientras el 
Jacior de umbral se ajuste de acuerdo a la media y la 
desviación estándar de la ventana vecindad sembrada. A 
medida que el tamano de la ventana se Incrementa. 
lexuras mas gruesas pueden sor omacas.en cuenta pera, 
expensas de reducción en certeza en lo que concieme 
la posición y delinción de las clases. La certeza a la 
membresía de una ciase se incrementa a expensas de la 
precisión espacial La exactitud en/a iasiicacióntodava 
está on consideración 

Tres aspectos son cricos en el crecimiento de 
regiones: el sembrado de un pixel. eltamaño de laventana, 
y el umbral. El prímero se realiza interactvamente al 
seleccionar un plrel representativo de la case de textura 
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a ser segmentada. Centrada en cada pixel sembrado se 
considera una ventana de tamaño 1 x 1 para la cual 
8 !/1! esuna primera aproximación de la dispersión de 
los valores de los pixels. Solamente ventanas impares son 
empleadas por lo que | = 2 + 1, donde a es un número 
patural. Los parametros a! y y) son la desviación 
estándar ylamecia dea población de pixels dela ventana. 
La localización de un pixel sembrado no es crucial 
mientras el umbral se ajuste de acuerdo a lo siguiente 





=,L .paralaciaso!, m 
n 


¡donde k,, es una constante que normaliza l umbral en el 
intervalo entre Oy 1. En otras palabras, para una ventana. 
ada, sil crecimiento de regiones es iniciado desde dos 
Pixels dilerentes pero que pertenacen ala misma clase de 
tono (Diferencia de Toro), se obtiene una segmentación 
basicamente equivalente siel umbral se apsta de acuerdo 
A las condiciones de la ventana local de acuerdo a (7) 

Es intervalo del umbral es entra O y 1. un valor céro 
"conduce al caso trval de una close vacia, el valor uno 
produce la imagen entera, la mejor estimación para uste 
Umbral será discutida después. Como'una primera 
Instancia, una estimación Óptima para el tamaño de la 
ventana se presenta para los tres predicados de 
untormidad empleados en osto trabajo. Construyamos. 
una distribución de frecuencias para los valores de los 
pixels (Gong y Howart, 1992) para cada ventana inicial de 
Clase de los pixole sembrados. Sea h(v) una distribución 
para la clase y para un tamaño de ventana 1 x 1, ondo 
ves el valor de un plrel y la función h(v misma os su 
recuencia de ocurrencia. Una gistancia entro un par de 
clases (1) puede ser ahora definida. 
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onde N, es el máximo posible valor de piel, en el 





prosanto baso, =255.La distancia promedio para todos 
los pares de clases os entonces 
NN, 
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donde N. es el número de clases, Para el predicado de 
uniformidad de Dierencia de Tono, el mejor tamaño de 
ventana 1 está dado para d máximo. 

Los predicados de unllormidad de Conelación y 
Entropía, especialmente adecuados para andisis de 
texturas, requieren de un enfoque diferente, SK) y (mun) 
son las coordenadas de los pixels sembrados para las. 
clases í y | respectivamente, la correlación entre las 
ventanas de las clases Z y 2 es 
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La correlación promedio de todos los pares de clases. 
esentonces 
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Elramano optimo de ventana es aquel que genera la 
minima Ó: 
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La entropla para la ventana de una clase dada es. 
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dona o, son los alamentos de la matiz do co-ocurrencia 
promediados sobre las ocho posibles, direcciones, tna 
cuenta detallada do esta fórmula sa da en a, siguien! 

occión. Para seguir el mismo razonamiento que arriba so. 
roquiero establecer un Valor de entropía para un par de 
ventanas de clasos. La matrz de co-ocurrancia so genera. 
sobra ¡el número normalizado: de: pares do, pixels que 
curan en una dirección, Para dos clases construyamos. 
una matriz, de co-ocurencia basada on los. pares que. 
ocurren phxela pixel on las dos ventanas, Para dos clasos: 
Idónticas, esta matriz do co-oourrencia encía solamente: 
un olomento diforonta de cero Igual a, uno, la entropía 
conjunta, sería coro; para clases diferentes, se obtandria 
'un Clorto valor alcanzando un máximo cuando las clases. 
soan totalmente diferentes. Por anto el valor promedio 
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1% máximo para ol mejor tamano e véntana 


Una estimación del umbral proviena de cada 
predicado de uniformidad y de candiciónos locales de las 
ventanas de los plxols sembrados, por tanto para una 
lorta ventana de clase con coordenadas centrales (, 
umbral para la correlación es 
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El umbral para la Dierencia de Tono está dado por la 
eouación(7) Las constantes, yl, sonpara propósitos de 
"normalización como ya se discutió arriba. El mejor umbral. 
se cotlene después de haber examinado la matrz de 
Jefes Malusita para todas ls distancias ntreclasos. La 
“ganeracion de un mapa de texturas requiero de una 
categorización global de texturas pará toda la imagen, 
stosahizoen dos pasos; una evalyación de entropía y un 
análisis de cúmulos sobre una imagen multespecral 
brida compuesta de dos bandas. la imagen original yla 
salda del analisis de entropía 
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5. SEGMENTACION DE TEXTURAS 
DE LA IMAGEN ENTERA 


Se desarrolló una segmentaciól nctependiente ce 
toxuras por medio dolcríteriode entropia(Haralck, 1970) 
basado en el concepto de la matíz de co-ocurrencia 
(González y Wintz, 1987), en combinación con el análisis 
do cúmulos, La maiz de coocurtencia se construyo 
ono sigue: Sea, Bun operador vectorial de posición, sua 
ke elnúmero permiido de valor de pixel enla imagen y sos 
kisuna matriz dx cuyo elemento (p.g) es el número de 
Vécos que un pixel con valor, ocurra en Una posición 
dada por Arolatvo a un piel convalor?,, Eloperador Bss 
¿define aqui para as ocho direcciones dadas por el Código 
de Frooman, Dado quo la correlación entre. pixels. 50 
desvanece rápidamento a mecácla que éstos e apartan, la 
longitud de Ano necesita sor más grando que la distancia 
nto vecinos cercanos. Entonces para un pel dado, sus. 
Vecinos con conectividad-8 son los que Usualmente 
«contlbuyen lo más ala mat. Notes que la cimensión 
de y, 6518 determinada por el máximo valor de pixel 
Indepondientemente del tamaño e, la imagen. En la 
prosante implantación, una vantana deslizante co amaño. 
vanablo (3X 3,5% 5... 11x 11) puede ser empleada para 
Calcular la matiz, asociada a cada vocindad de pixel. El 
elemento m,, de y almacena ol número do veces que un 
par de valoras (4) ocurren en la posición rolativa dada: 
por una voz evaluada y,, Cada blemontomm,, os divido. 
por el número lotal de pares, el resultado 68 una nueva. 
atra Cllamada matiz do co-ocurtonclo, que Garactoriza 
on buena aproximación a probabilidad conjunta de que un 
par de pirols que sals/agan tengan valores (+). En 
lras palabras la matriz Cs una aproximación la función. 
de densidad conjunta. Exito simstra axial un las ocho 
direcciones escogidas de lal manera que Cs una matriz 
simétrica. Esta matriz.as de, hecho el histograma 
bidimensional cuando la ventana Iguala al tamaño de la 
imagen. 

Para caraciorizar el contenido de la matriz co 0o- 
ocurrencia, se ulza l crio de entropía 
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donde, sonlos elementos do la matiz C, Esta expresión 
está basada on loyos tormodinámicas y proporciona una. 
medida de orden y probabilidad de un sistema cerrado de 
la siguiente manera. 


Entropia Ata - Menos Ordén - Ata Probabilidad do Ocurrancía 
Entropia Baja - Mayor Orden - Baja Probabilidad de Ocurenci 





Lalunción de densidad conjunta c,, es siempre menor 
o igualauno. Cuandolos valores delos phkels.onlaventana, 
deslizante son todos Iguales, una de las c,,= 1 y el esto 
on cero, de aquí que Entropia = O. su minimo valor. Silos 
valores de los pixels son todos diferentes, entonces cy < 
1yEntropla>0, conunciertomáximovalor.Por tanto Una. 
área de textura suava falto orden) en la imagen le 
corresponde un valor bajo de plrel (tano oscuro) en la 
imagen detextura, Una textura gruesa se traduce en tonos 
úclaros,estoimplica que un mapa de texturas segmentadas 
es posible al formar cúmulos de tonos (valores de pixels) 
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dea imagen de textura. Esto so realizó por medio de la 
siguiente hipótesis: la Imagen original contiene la 
información espectral (tono) de la escena y la imagen de 
textura leva la, información espacial, por tanto una 
combinación de estas los imágenes en una imagen 
multespectral de dos bandas debe proporcionar al 
analisis de texturas de la escena. El análisis fue levado a 
cabo con al Algoritmo de Cumulos amado de isodata. 


6. RESULTADOS 


Dos imágenes SAR de 512 x 512 pixels con una. 
resolución de 25 m fueran ullizadas para probar los 
métodos desarrollados en este trabajo. Una es una 
subregión de una imagen SEASAT conocida como del 
anticinal de Obayos, localizado en el norte de México, la 
tra es Una imagen serotransportada del norte de Ho- 
lana, Fue con!imado expermentalmente que la función 
de densidad del "specke” para la imagen SEASAT es 
Gausslana y Rayleigh para la aerotransportada, A ambas. 
imágenes se aplicaron una serie de fitros (sección 3) con 
losrosultados arriba ciscutidos. Aun cuando elfiro de Lee 
muliplicativo borra algo los bordes de la imagen. su 
comportamiento. tan bueno como+l del ftro geométrico 
on lo que al cracímianto de regionos se reñiere. Para los 
propósitos de clasificación de texturas basada, en la 
entropía, os mejor ullizar el fitro geométrico de antemano, 

Se aplicaron dos tipos de procedimientos de 
crecimiento de regiones a ambas imagenes. Para la 
imagen SEASAT, sa distinguleron claramente dos clases 
de texturas, Wtlizando. al criterio de la correlación 
hormalzada y con una ventana de 5 x 5, estas clases 
prosumiblemente corresponden a dos unidades 
Iológicas, se requiera de más trabajo en esta dirección. 
Para la Ímagon aerolransporiada (Fig. 5) so sapararon 
claramente una cludad y un. campo de cutivo empleando 
bl mismocrierio yla misma ventana. Algunas pruebas con 
el criterio de entropía fueron también ejecutadas y una 
vantana de 5 x 5 pareca proporcionar una segmentación 
más detallada para la textura, aún quando el costo 
computacional es mayor. Todas las pruebas fueron 
aplicadas dospués de una iteración del litro geométrica. 
Sa trabajó una ciasicación separada por modo de la 
matiz de co-ocurrencia y el criterio de entropía para toda: 
la Imagen SEASAT. Una operación subsecuente de 
cúmulos con 4clases indica una clara categorización (Fg 
6) de texturas. Las baridas empleadas en esta operación 
de cúmulos fueron la imagen filtrada con el filtro 
¡peométrico con una eración (a= 10%) yla que rosuta de: 
la medida de entropía. Se realizó una comparación con 
una operación de cúmulos (Fig. 7) utizando la imagen 
ciginal sin frale y la que resulta de la medida de la 
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“antropia. En ambos casos se evaluaron. la matriz de 
covarianza y la distancia de Jefires-Matusa (Tablas 1y 
2) Una inspección visual (Figs. 6 y 7) y la comparación de 
estas matrices muestran cjaramente una mejoría en la 
clasificación cuando la imagen fitrada es incluida en el 
análisis de cúmulos. 


7. CONCLUSIONES 


Se ha presentado un método para detectar objetos 
texurales especificos o para generar una segmentación 
¿lcbal en imágenes SAR. Se ha proporcionado tambien 
una metodología para la reducción del "9peckle" y su 
Impacto enel analisis do texturas. Se presentan resultados 
empleando imágenes SAR aerotransportadas y 
satelitales, pero elmétodo es aplicable a cualquier Imagen 
de rango generada por medio de radiación coherente. 

Encomparación con trabajos previos (Connets.l a, 
1984; Pecie y Franikin, 1991; He y Wang, 1990), elcáiculo 
dela entropía aquíestá basado en el promedio de la matriz 
de co-ocurrencia promeciada sobre las ocho direcciones. 
definidas por el Codigo de Freeman, Esto reduco la 
"dependencia de la textura sobra posibles variaciones en 
condiciones de luminosidad u otros corrimientos 
monctónicos de los valoros de los pixels an la Imagen 

La detección de objetos. loxturalos. uilizando el 
procedimiento de crecimiento de regiones puedo ser visto 
Como un paso previa para evaluar caracioríaticas 
morfológicas y Jopológicas de aspectos selectos de la 
escena. Estos datos pueden sor, a su vez, empleados 
como entrada a una base de datos para el manejo do 
recursos relacionados con la escena. 

Laroducción del"specklo" llano claras consecuencias 
sobre la segmentación de texturas; existe Una mejoria 
visual después del fitaje, resultando en un mapa suavo' 
(Fig, 6) de texturas, Antes del Majo dol "gpecklo" se 
observa una fragmentación (Fig. 7) de las clases do 
texturas. Además de esto, la matiz de covarianza y la 
distancia de Jefínes:Matusita muestran una mejoria on 
cuanto a la compactividad y separación do clases. 

Una generalización da la evaluación de la ontropía a 
Imágenes muliespectrales podría proporcionar un 
clasificador contextuakespectral, no obstante los cálculos. 
serían más bien elaborados involucrando un alto costo 
computacional 
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TABLA 1. MATRIZ DE COVARIANZA Y DISTANCIA JEFFRIES - MATUSITA ANTES DE LA REDUCCION DEL "SPECKLE"_ 








Clase de Textura Matriz de Covarianza Distancia Jeffries - Matusita 
Suave (1) 10025 32755 Y 2 3 4 
32755 — 151532 
Regular (2) 121822 200924 a es A 
2000.24 363645 2 0779 0000 1989 1.999 
Intermedia (3) — 25804 -2733 > Ea 17 00 a: 
273 09216 4 2000 1999 1563 0000 
Rugosa (4) 247.10 15548 
15548 — 78823 








TABLA 2- MATRIZ DE COVARIANZA Y DISTANCIA JEFFRIES - MATUSITA DESPUES DE LA REDUCCION DEL "SPECKLE” 














Clase de Textura Matriz de Covarianza. Distancia Jeffries - Matusita_ 
Sa AA 
e e 1 0.000 1.015 2.000 2.000 
roga seno ar 
IEC 
A 
E 
hugo mas tasa 
dE 
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Abstract 


Data integration ls obviously nota new concept and has been pursued for many years 
¡on an analog basis in many Earth science disciplines, However, rapid advancos in image 
analysis hardware and software have allowed for greater flexibllty and innovative toch- 
hiquos for combining and integrating digital data. 

Many techniques exlatfor combining digital data. However, the end products (Lo., color 
“composite images) are ofton difícult to interpret quantitatively and qualitatively as tho 
statistical properties of the data have been manipulated and, thus, tho original Integrity of 
e data ls not left intact 

“The color display transtorms such as intensity-huo-saturation (HS) can bo usod to 
produce more offectivo and controlled visual presentation ofthe data forboth qualitative and 
ualltativo Interpretation procedures. Tho IHS color transtorm has scon many applications 
1or the display ol remolely sensod data, 

"This papor describes examples of INS transformed images for tho Integration o TM 





LANDSAT and SAR - ERS=1 data. 








1.- INTRODUCCION 
nadelas mayores imitaciones que so encuen- 
tran en el uso de la Percepción Remota os la 
presencia de nubes, En un momento dotermi- 
ado una gran parto de la superficio torreta 
stá nublada y on algunas rogionas talos como los tropk 
lugar donde so encuentra situado nuostro palo, la 
eoberlura es casi pormananto, 

En consecuencia, dobido a las conciciones atínostó: 
/icas y a la nubosidad, un enfoque mulitemporal para al 
mapeo del uso de la Narra sulta de ciertas limitaciones, 
Esta problema ha sido parolalmonto resuelto por ol Sisto- 
¡ma SPOT y su capacidad de observación fuera del Nadir 
y por el uso combinado de imágones SPOT y LANDSAT 
paralamismaregión, pero'n algunas áreas tropicalós aún 
hay una elevada carencia de información y es all dono la 
Imagen aciar y au potencial para el mapeotomático pueda. 
ser de especial importancia. 

Porotrolado,las técnicas de procesamiento digita dé 
Imágenes han evolucionado an elsentido da integrar datos. 
de diferentes sistemas sensores 

Bajo esto contexto, este artículo eno ol propósito de 
presentar un proctueto, resultante de la integración de una. 
Imagen radar SAR - ERS-1 y úna Imagon LANDSAT TM. 


2.- ZONA DE ESTUDIO 

















Luego dol andísis de la disponiblidad de datos 
orbitales de los dos sensores, tomando en cuenta que la 
información radar se encuentra menos distorsionada en 
zonas planas, se seleccionó la zona de manglaras al Sur 
de la ciudad de Guayaqui, y el Valla de fos Chill, 
considerando además la importancia ecológica y econó» 
mica de estos sectores. 





« PRODUCTOS SATELITARIOS UTILIZADOS 


TMLANDSAT 
Imágenes gocreléronciadas, adquiridas el 15.40 ocu 
bro de 1991 (Imagen Quito) y 29 desabrl de 1991 (Imagen 
Guayaquil), bandas; 3(0.63-0.69 um). 4 (076 0.90 um) 
y5(1.55-1,75 um): Rosolución espacial 30 mts. x 30 mis 


SAR -ERSA 

Imágenes adquiidas, ol 17 de abel de 1992 (Orota 
13949 Guayaqui) y al 19 do mayo da 1992 (Orbita 4407 
Cotopaxi) banda G (longitude onda de 5,6 cm, aquiva- 
lente aunalrecuenciado5.3GHZ ),resoluciónespacialdo. 
30m. x30mn..AnigulodoIncidonciacdo23*, polarización VW, 
contitraje basado sobre la mediana, el cual se haaplicado 
'conéxitoparareducirel moteado brllane as imágiones. 
racár (Blom y Dal, 1982) debido a que suprime eluido. 
respetando los bordes. 


SOFTWARE 
MERIDIAN versión 312, 


4.- METODOLOGIA 


La metodología empleada consista de dos fases; 

REGISTRO GEOMÉTRICO delos datos delos produc 
tos a serintegrados, FUSION DE LOS DATOS registrados. 
'en un producto único, 





« REGISTRO GEOMETRICO 


“Tanto a imagen TM LANDSAT como la imagen SAR 
ERS-1, fueron corregidas geométricamente con canogra- 
a básica a la que se pretendió ajustar las imágenes 
separadamente y luego ajustarias entre sí, considerando. 
¡como referencia la imagen TM LANDSAT. 
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Para calcular cualquier lunción ente 008 0 más varía 
blas, es necesario contar previamente con una seno de 
obsorvaciones comunes a todas ellas. La calidad del 
ajuste dependerá dela prectsióncon que selocalcenesos 
puntos, el establecimiento de puntos de control resulta la 
lase más crucial del proceso de corrección y la que 
demanda mayor dedicación humana. 

Para que el ajusie entre imagen y mapa. y entra 
Imágenes. sea conecto, se requirió atunder a tres aspec- 
1os en la selección de los puntos de conta! número, 
localización y distribución. 

El número idóneo de puntos de control depende del 
Iamañoy de latopogralia delterreno.Silas deformaciones 
que presentan zon importantes, conviene aplicar funcio 
nas de taraformación complejas. Cuanto mayor sea al 
grado de la ecuación de ajuste, se precisa tambien un 
mayor número de puntos e control 

'Enlas zonas de estudio el contrasto topográfico no es 
do 1YaKGondongia, sin ambargo se seleccionó un polinomio 
do translormación de negundo grado en la segusento forma: 





Fay) 0 AAA ae ay + ay 
y aproximadamente 30 puntos de control para cada regi: 
I1o(matomaticamanto son necesarios solamente seis pun» 
109) para garantizar una buena gontormicad. 

En cunIvo a la localización, sw deterrinaron puntos. 
ú'latamento Igentlicabies, profontiemorto rasgos numa, 
hor no sujetos a dinamismo Temporal (piscinas ce las 
tamaroneras. corta do caminos, puentes, etc.) 

Rospacio a la dixtibución, las puntos de conto! sa. 
ústuaron uniformemente sobro lalmagon Esto ovio orrores. 
debido a una ponderación uxceniva ce algún aoctor del 
Aroa de estudio. 

La bondad cel grado de ajusla se mide con los 
rosiduales, Como 04 Dion sabido. el residual do la rogro- 
ion ex la foronala entre el valor estimado y el valor roms, 
para coda una do la» Obuervaciones empleadas on el 
proceso, En nuestro caso se aptó por un residual maximo. 
do 1 pino! (30m. 

En primera Instancia el residual era supero al inicia: 
manto marcado enaltrabajo, se procedió suprimir elolos 
puntos que oleciar una mayor desviación on las estima- 
lonas, Una voz olíminacios, so calcularon de nuevo los 
coofIciontes 

Lon polinomios de taraformación ctados anteror- 
mante perrton vaduci coordenadas de la imagen TM 
LANDSAT a la imagen SAR - ERS -1. Ahora Dion, sí 
Fretonde crearse una imagen que se corresponda ade: 
cuadamente son esas coordonadas, resultaba necesario 
lrasvasar de alguna forma los niveles digitales 00 la ma- 
¡gen original a esa nueva posición. Para este electo se usó 
la nlerpolación Diineal, la cual promedia el nivel cigital de 
los cual plkoles más cercanos en la imagen inicial 








'USION DE LOS DATOS 
REGISTRADOS: 


TRANSFORMACIÓN HS (Intansiy. hue, saturation) 

El sistema de color ROJO-VERDE-AZUL (RGB por sus 
'nombres en ingles) comunmente usado para desplegar 
ves canales de las imágenes sateltanias en Percepción 
Remota, es simple y efectivo, sin embargo tiene algunos 
electos, entre los que se destaca que el sistema RG8 no 
úestábasadoan atributos del color definidos tácimente yen 
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consecuencia las variaciones del color dolinidas por la 
'mezcia de las colores primarios no pueden ser descritos 
Tuméricamente, 

Un sistema efectivo que supera todas ostos inconwe- 
úentes es la ranstormación IHS, la cual esta delinida por 
11es atributos de color, Independientes, ortogonales y 
lácimente percibidos: Intensidad, Matiz y Saturación 

El sistema de coordenadas IHS puede representarse 
'geoménricamente como una esfera, 

La intensidad, la cual representa la energia lotal o la 
ú'rllantez de la imagen, está deliida por el radio de la 
sera. El matiz representa la longitud de onda pramecio 
¡del color y está definido por el ángulo centra dela esfora, 
La saturación 05 considerada como la pureza del color 
(porcentaje o lut blanca enla imagen) y astá detinica por 
la colaLIUO do la estora 





7.- FUSION 


En el presente trabajo, las bandas TMS, TMA y TM3 
lueron procesadas por latranstormación del ospacioRGB, 
al espacio 1HS. 

El histograma de la componente $ (saturación), luo 
¡modificado con el objetivo de adicionar valoros para tomar 
los colores más vivos en el producto final, 

La componente | (inmensidad) fue susttuida por la 
banda G del ERS-1 

Finalmente tue realizada la transtormación IS inver 
sa. obleniéndose Como imagen resultante una compost: 
ción falso color que reúna las caracteristicas da LANOSAT 
yERSa 

La guiente figura lustrael proceso detranstormación 
para la combinación delas dos Imagenes, 


ANALISIS DE RESULTADOS 


'Unabuena precisiónenel ajuste fuel objetiva de esta 
investgación, con la metodología anteriormente clíada se 
halogrado en buena parte esta meta. El numero de puntos. 
e contras la precisión con la que fueron determinados. el 
grado del polinomio de ranstu. mación yla suava topogra- 
fía contribuyeron para que setenga un buen fegistrocamo 
producto final 

La magen falso color, obtenida con la integración de 
las bandas de la imagen TM LANDSAT y la banda C de la 
imagen SAR - ERS-1 reune las características de las dos 
imágenes componentes. La gran ventaja de esta imagen 
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Imagen oríginal TM, 
composición bandas 4,5,3 (RGB), 
zona de Guayas. 







'Sinergísmo LANDSAT TM/ERS-1-SAR 
transformación (RGB) / (SAR 1.H.S.), 
zona de Guayas. 
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Sinergísmo LANDSAT TM/ERS-1-SAR 
transformación (RGB) /(SAR IWS). 
sector de zona de Guayas. 





Sinecgísmo LANDSAT TM/ERS-1- SAR 
E transformación (RGB) /(SAR HS), 
zona de Sangolquí. 
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es que simultáneamente elía ofrece. por una pare, las 
propiedades importantes de las microondas, tales como: 
Penetración en las nubes, penetración en la lua incde- 
pendencia de energia emitida por elso!, penetración enlas 
Copas delos arboles y una interacción con la matera en 
Junción de las propiedades del terreno (constante 
dieléciica y rugosidad de la supericia) diferente a la 
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interacción de las longitudes de onda del visible y del 

intrarroo usado en lasécnicas tradicionales de la percep- 

¡ción remota. Y. por otra parte, el poder espectral de las 
“TM LANOSAT. 

'El producto resultante del sinergisho TM LANDSAT/ 
'SAR ERS-1 ademas de la faciidag que permite extraer 
¡oda información de un único producto, es mejor para la 
Interpretación visual que se realiza de manera más preci- 
sa, en este caso, en lo que se rellere a la infraestructura 
'camaronera, lo cual permitrá una Identificación exacta de 
les piscinas. 

La información SAR ERS-1 coadyuva para este fin 
debido a que las paredes de las piscinas se comportan 
'comorefieciores de esquina, es decir, dan lugar areflexio- 
nes especulares en dirección a la antena, originando 
puntos brilantes que permilen ta fáci identlicación, sin 
recta ls concisones de rubendas, 

La respuesta espectral de LANDSAT sumada a la 
úhabiidad del Radar de Apertura Sintética para discriminar 
entre cuerpos húmedos, agua y superficies secas hace 
¡que este producto en el futuro se convierta en un instru. 
mento adecuado para el uso de estos datos en el estudio 
¡e las condiciones de humedad del suela. 

El producto final guarda las características 
espectrales de LANOSATe incrementa la caracterísica de 
rugosidad propia de RADAR, resuelve el problema de la 
nubosidad. en la zona urbana permie distinguir las áreas 
de mayor densidad de construcción, permite también 
'discrmnarlascerreteras. las cuales se encuentran realza 
das. 

Los lineamientos, geológicamente hablando, se pre- 
úsentan con mayor clardad en la imagen resultante, pro- 
¡acto de la aportación del relleve de los datos SAR ERS- 
1 y de la caracteristica del color de la vegetación de los 
tatos TM LANDSAT, 

Gracias a la resolución y a las caracierísicas de 
terodispersión de RADAR, es fácil Idamiíicar los barcos 
¡que se encuentran en la zona de estudio, como se puedo 
observar en la siguiente fotografía. 

En lo que tene relación con la zona de estudio de 
Sangolqu no se obtuvo una buena superposición de las 
¡rmágenes debido a la Iopografía del terreno, que a pesar 
de ser un valle está distorsionado por las elevaciones 
úcscundantes.nlairmagenRADAR. Snembargo se obluvo 
úl siguente producto en el cual se puede visualizar el uso 
del suelo en zonas de fuerte refeve: 

Actualmente, en varios centros de Investigación, se 
están realizando trabajos para lograr un mejor procesa- 
miento de las RADAR con la aplicación de 


imágenes. 
sofisticados fitros para eliminar el SPEACKLE y el ruido y 
con modelos de terreno para mejorar la corrección 
¡geométrica en relieves pronunciados. 

Se espera que estos tipos de aportes coadyuven para 
el mej aprovechamiento de esta Iécnica de registro en 
zonas de fuerte topografia 
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ERS-1 : EL ESPACIO AL SERVICIO 
DE NUESTRO AMBIENTE, 





Laurent POLIDOR! 





AEROSPATIALE, Cannes (FRANCE) 


cludad de Cannes (Francia) el primer 

simposio dedicado al satélte ERS-1. Para 
cerca de 400 Investigacores, principalmente europeos. 
esto evento fue la oportunidad de Intercambiar ideas y 
resultados despues de un año de funcionamiento del 
satolte, Sacándose conclusiones sobra la calidad de los 
datos proporcionados hasta ahora por los diferentes 
¡sensores dol ERS-1 y más que todo por el SAR (Radar de 
Apertura Sintética). Por otra parto, resultaios sumamente 
prometedores Iueron presentados en varios campos de 
aplicación de estos datos, Estos resultados provienen de 
studios experimentales lavados a cabo por laboratorios 
úBuropeos de Investigación y patrocinados por la ESA 


CALIDAD DE LOS DATOS 


Según las conclusiones de la sesión dedicada a la 
calibración del SAR, la calidad de las imágenes, tanto 
radlométicas como, geomáticas, es casi perfecta. Los 
roprosontantes deja ESA presentarontos resultados delas 
mediciones lovadas a cabo durante un año para estudiar 
la evolución del nivel de ruido y la estabilidad de los 
ilerantos parámetros instrumentales. Los usuarios tam- 
bién llenen oxigencias en cuanto a la calibración de los 
datos, y los investigadores de la ESA y de los laboratorios 
asociados intercambiaron informaciones. sobre la 
calibración relativa (el SAR debe dar la mismarespuesta 8 
Observa la misma superficio) y la calbración absoluta (el 
valor de esta respuesta debe coresponder exactamenta: 
a la energía rellejada por la superficie), Se usaron 
instrumontos de calibración, pasivos y activos, para los. 
Cuales se conoce exactamente la respuesta eórca, así 
que se pueden hacer correcciones en caso de pequeñas 
dleroncias. 

Al finalizar la sesión se do en broma que gracias al 
ERS=1 hablan calibrado los instrumentos de calibración. 


ol mes de noviembre de 1992 la ESA 
(Agencia Espacial Europea) organizó en a 


APLICACIONES TEMATICAS DEL SAR 


Oceanografía 


El satélite ERS-1 fue diseñado, más que todo, para. 
aplicaciones oceanográficas esencialmente, y muchos 
resultados fueron presentados en este campo. Modelos 
anteriormente desarrollados y validados con el SEASAT 
fueron refinados con las Imágenes del SAR para 
Investigaciones de oceanografía fisica relacionadas con 
tonómenos muy variados como corrientes, remolinos, olas 
Internas, vientos de superficie. Varios laboratorios 


Interesados en modelos climáticos globales. pudieron 
sacar provecho de la sensiblidad del SAR al estado de la 
superficie del mer. 

Aunque las ondas de radar no penetran en el agua, 
permiten detectar bancos de arena en el londo del mar por 
¡el cambio que producen en las características do las olas. 
El fenómeno ya se había observado con las imágenes del 
'SEASAT. Las imágenes permitieron llevara cabo estudios 
de batimetria y aplicarlos al estudio de fenómenos 
costeros (geomortología, hidrología). 

Fuera delas investigaciones puramente cientificas, los 
datos del ERS-1 fueron empleados para aplicacionos 
operacionales tales como el. contro! de la contaminación 
(manchas de aceite) o del tráfico de barcos (el SAR es muy 
sensible a la estela de los barcos), En la figura 1, donde. 
cada color (azul, verde y rojo) representa una lecha 
lerente, la superposición mulitemporal de 3 imágenes. 
¿el SAR permte distinguir limites de predios agrícolas. 
variación de mareas y otros fenómenos. En la figura 2, 
correspondiente a una sola facha con colores reales. 
simulados, se aprecian detalles agricolas, urbanos y de: 
¡oleaje.en arma muy niida, Ambos ejemplos corresponden 
al Norte de Holanda, 


Cartografía y Geociencias 


Debido a sus caracteristicas instrumentales y sobro 
todo a su ángulo de incidencia muy estrecho (23") el SAR. 
del satólto ERS-1 os muy sensible lroliovo terrestro, Esta 
sensibilidad consiste en fenómenos de deformación que 
pueden ser considerados como defectos por los usuarios 
interesados en el estacko dela suparfcio ya que vuelven al 
análisis más dificil, Pero son muy valiosos para aquellos 
usuarios interesados en observar, medir y mapoar las 
formas del terreno: ingenieros civiles, geólogos, 
cartógralos. 

Algunas aplicaciones geológicas de las imágenes del 
SAR fueron presentadas en el simposio: estudios 
gscmorfológicos y mortotectónicos en zonas de selva 
tropical (Guyana, Venezuela); análisis de procesos 
volcánicos en islas del Pacifico: observación del estado de 
superficies en zonas desérticas. 

Los resultados de la llamada técnica de in- 
terteromenía, tanto para aplicaciones topográficas como, 
geofísicas, tuvieron mucho éxito. En efecto, la 
comparación de imágenes del SAR tomadas desde 
órbitas muy cercanas permiten detectar variaciones 
espaciales o temporales del orden de la longitud de onda 
($ em). Fuera de la capacidad de construir modelos 
¡digitales de terreno con una precisión de algunos metros, 
so tiene como consecuencia la posibilidad de detectar 
movimientos centimétricos afectando áreas extensas: Las: 
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aplicaciones ambientales (observación de desliza- 
mientos) y tectónicas (medición de movimientos de 
placas), pero también volcánicas, hidrológicas y 
¡glaciológicas, apenas se están vistumbrando, y muchas 
investigaciones quedan por desarrollar 

Estudios Forestales 


Con sus parámetros instrumentales únicos 
(Nncidencia, frecuencia, polarización) el SAR del ERS-1 no 
Tuo diseñado para estudiar la vegetación, pero permitó 
extraer bastante Información en este campo también. 





para evaluar la adecuación del SAR para estadisticas 
Agricolas y predicciones de cosecha. 

Sin embargo, la mayoría de las aplicaciones rela- 
elonadasconla vegetación son forestales. EISAR fue usado 
para la detacción de cambios estacionales en zonas de 
selva y para identificar diferentes ciases de árboles. 

Entro los estudios forestales, el más ambicioso a 
Interesante es bl proyecto TREES (Tropical Ecosystem 
Environment Observations by Satelites). Su meta es la 
dentlicación y clasificación de especias de árboles entas 
¡selvas tropicales. El nterés de este proyecto es que todos 
los aspectos que puedan inf en la confiabiigad de los 
resultados. son tenidos en cuenta. como los procesa- 
mientos raciométricos y geométricos anteriores al andisis, 
y el nivel de proprocesamiento que sa debe exi 


LOS OTROS SENSORES 


Acdomás col Radar de Apertura Stética, al satalto 
ERS-1 está lovandouna carga úl compuesta de un radar 
alÍmetr, un radar disporsómetro y dos raciómatos pa- 
vos (uno lórmico, el tro de microondas), Como en el 
080 del satólto SEASAT, todos estos sensores fueran 
diseñados en primer lugar para apícaciones oceano. 
¿rálicas, poro proporcionaron también informaciones muy 
Valosas para aplicaciones continentales. 


El Radar Altímetro 


El radar alimetro mide la altra del satelte con 
(ospectoala superficie del mar. Conociendo la rtta, sta 
medición de alura permite determinar la topografía del 
mor con una precisión del orden de 10 cm (en promedio 
sobre droas de varios klómetos de extensión), Como se 
hablahecho con elalimetro de SEASAT, se sacó provecho 
dolasenslblidad del adaralaaltura media (estudio delas. 
mareas), a la pendiente (corientes, remoinos y oros 
fenómenos de circulación) y ala rugosidad del mar (olas. 
viento...) Ealímeto es muy adecuado para el estuco de 
onómenos de poca amplitud pero que se manifestan 
sobre áreas muy grandes. Sus datos permiten almentas y 
validar modelos cimatológicos a escala planetaria. 

Aunque fue desarrollado para estudios oceano 
grálicos, el radar alímatro fe utlizado para aplicaciones 
topográficas también. Su resolución espacial no permite 
resilul las gooformas pequeñas, pero las frecuencias 
bajas de roleve son representadas con un alto grado de 
confiabilidad. Ha sido posible realizar modelos igiales 
del Continente Antático calculado con el radar alimetr. 
Somejantos mapas pueden ser uilizados para com- 
plomentar le cartografía de regiones donde sigue 
haciendo ala, o pararegisrar datos de mayor detalle pero. 
ral localizados 
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El Radar Dispersómetro 


¡Con sus 3 antenas, el dispersómeto fue ciseñado para. 
observar las olas desde 3 drecciones complementarias de 
manera casi simultánea. Anteriormente, otros disp 
sémetos, en paricular el de SEASAT, hablan demostrado su 
“aptitud para medi la velocidad y la iección del viento en la 
superfice del mar, Por lo tanto, varios laboratorios de 
meteorología están involucrados. en Investigaciones 
“aplicando los datos del dispersómetro del ERS-1, Resultados 
¡prometedores fueron presentados en el simposio. 


Los Radiómetros Pasivos 


Además delos sensores activos, elsatélte ERS-1 tiene 
'dosradiómetros pasivos. Uno funcionacon microondas, el 
¡oo en el infrarojo térmico, Aparte de su uso para las 
¡Correcciones atmosféricas de los datos del allímetro, los 
radiómetros tenen varias aplicaciones temáticas basadas 
'enlas mediciones de temperatura y humedad en la super 
ficiedel mar o dela era. Varios resultados fueron presen- 
tacos sobre la validación de modelos físicos, usando los. 
datos de los raciómetros del ERS-1 porn lado y medicio: 
nes de campo o radiosondeos atmosféricos por otro. La 
figura 3 muestra algunas aplicaciones de modelos de 
temperatura superficial, 


PERSPECTIVAS LATINOAMERICANAS 


La mayoría de las aplicaciones temáticas 
consideradas en este simposio pueden contribuir a los. 
estudios ambientales en América Latina. En efecto, las 
:igaciones llevadas a cabo por laboratorios 
úcienticos están demostrando que el ERS-+ ha dejado de. 
ser estrictamente experimental para volverse una 
herramienta operacional que puede ser utlizada en casos. 
reales. Además, nuevas aplicaciont están 
vislumbrando con la posibilidad de medir parámelros. 
¡eofísicos y de hacer mediciones de desplazamientos a 
úivel de la longitud de onda. 
Paraunusooperacional,otros factores, en particul 
trecuencia de las tomas y la disponibilidad de 
estaciones receptoras, tienen que ser considerados, y 
todos los aspectos del sistema (sensor, órbita, recepción 
y preprocesamiento, precio de los datos) deben ser 
analizados porlos usuarios parala programación eficiente 
¡de una aplicación operacional de los datos del ERS-1, La 
figura 4 muestra un sector de imagen adquirida desd la 
Estación Alemana en la Base Antártica Chilena O'Higgins. 
lo que posibitará importantes estudios an el Extrerno Sur 
de Sudamérica yla Antártica 


CONCLUSION 


Aparte del SEASAT, que no duró más de 3 meses enel 
año 1978, el satólte ERS-1 está llevando el primer radar 
espacial. Cuando fue lanzado en julo del 1991, la 
selección de una nueva longitud de onda y el proyecto de 
usar los datos a nivel operacional constluian un reto 
grande. Ahora, toda la comunidad cientifica está 
¡Convencida del interés del SAR y de los otros sensores 
para aplicaciones operacionales tan variadas como los 
Inventarios forestales, elmapeo de losriesgos nalurales, el 
Control del tráfico y de la contaminación en el mer... La 
¡conclusión de este simposio fue que el ERS-1 esrealmente 
un satélite al servicio de nuestro ambiente, 
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DIMENSION...el receptor GPS compacto con 
precision igual a las dimensiones de un centavo 
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Intero provides a wide range ol remote sensing technologies trom 
data collectlon throwgh analysis, including; traíning, leennotogy transfer 
end consulting services. 


Spociollles include: 
Imaging radar tor 
Geological and Forestry Appilcations 
Land Use/Agriculture Mapping 
Cartographie mapping 
Training programs 
Rader Simulator for all spaceborne tadar plallorms including ERS-1 and Redarsar 
Forest Penetrating Radar lor prollles beneath dense forest canopy 
Environmental Surveillance (Oi! Spills, Flood) 
Environmental Impact Studies 
Westher monitoring and modification 
Geographic Information Systems 
Image Processing syatems 
Systoms Englnocring 


els an edo 
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The Systems Engineering Company 


Systems for: 
Earth Observation, Environmental Monitoring, 
Resource Management and 
Geographic Information Management 
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ThePerfect Partnership 


Of High-Performance Imaging 





The future belongs to those who define it 
and the future of graphics and imaging systems 
has been defined by Intergraph* and VITec 
Announcing the Intergraph ImageStation 6187”,a 
new standard for RISCbased workstations that 
tightly integrates the worlds of sophisticated graphics 
and high-performance imaging, 
ImageStation 6187 
The ImageStation 6187 integrates ViTec's imaging 
chipset with Intergrapt's EDGE Il 
graphics subsystem establishinga 
new performance benchmark for 
graphicsand imaging workstations. 
For Application Developers 
andusersin CIS, C3l, Medical and 











AndSophisticated Graphics. 


INTERGRAPH 
€ 
VITec 


The Imaging Alliance. 





Scientific Visualizatlon markets this means unprec 
edented processing speeds ol 25,000 polygons per 
second, 300 MOPS (Million Operations Per Second) 
and18 MFLOPSon rasterimagedata. 


Open Environment 
The Imagestation 6187 is fully compatible with all 
other Intergraph RISCsystems,and incorporates 
available industry standards such asthe X Window 
System”, the UNIX” Operating System, and leading 
networking protocols. 
Formoreinformation-ora 
lookat thefutureo! workstation 
technology-call Intergraph 
at8008263515017722700in 
Huntsville, Alabama. 


AG 0 





PATROCINIO IMEROAMERICANO E. INTERNACIONAL 


SELPE! 


VPSTA software de percepción remota 


y modelamiento espacial 


VISTA! es un con yde módulos de software de 

proce ma, creado. 

para satisfacer las ne los protesonales en manejo 

¿le recurs tatuados, gcologísy control ambiental 

¿Con 20 años de experiencia desarrollando y manteniendo 
andare y sistemas completos de perecpeión 

OS disonó VISTA especificamente para logra l 
























ruáximo de rendimiento en evtaciones de trabajo y 
computadores de somsolahasados on UNIX, 
X-Windows 

Silicon Ciraphica*, MY Apollo" y 








e gráfica para ol usuario po apurte/Mlspare 

ja con reparadores de bandas de 8624 bits 

+ Tamaño de imajen nicas lación (00m) y selección 
¿le Aca on) vial 

+ Anotación/edición iractiva de información raserivector 


encamiick 














Nperpueata en imúgenes raster 
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+ Incluye formatos comunes de intercambio de daros 

Compatible con un amplio rango de dispositivos de 
entrada y salia (1/0) 

+ Sobre 300 funciones algoriunos) disponibles, incluyendo 
calce de imágenes. funciones geométricas, clasificación. 
modelsmiemto espacial, Altos, manejo de información de 
vectores, apoyo SIG y presentación de datos 

= Acepta información satelital en cualquiera de lor 26 formatos 

= Mantiene 34 tipos de proyecciones de mapas 

+ Móxulo de modelamiemto raster SIG. 

= Completo apoyo al usuario y un año de garantía incluido 
enel precio de venta 

















Para mayores detalles, lame escriba a: Itermaional 
Imaging Syaem, 1500 Bucheye Drive, Milpias, Calforía 
195035. Tel. 4ORJ412-1400, Fur ADR/AJ3.0065 


+Marcas Registradas de sus respectivas compañías, 
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Establecimiento Alemán de Investigación Aeroespacial (DLR) 








Arva Espacial 


Transportadores Espaciales 

+ Vuelo Hiparsórico y Problemas 

d Foonitais en la Almónlora 
Oolalos 





Apicación do Tacnclogías Espaciales 
+ Comunicación y Téenican 
de Localación 
+ Sislomas de Exploración 
«Exploración dol Mociowmionto 
« Clenca do los Matailos 
y Clancias Bilégicas 
Totana relacionada 
con Mislonon Espacialos 
' Coniron Espaciales 
 Sstomas do Conti on Tora 
* royoción de Mntonos Espúcialos 
Tareas organizalwns 
ono Hon osacia 
“a Eocución do Proyoctos vinculados 
¡on Cionclas y Tacnicas Espacial 
+ Ovas areas cIganizatwos 














Area Asrondutica 


+ Navegación y Seguídad do Vuolo 
+ Tecnclogía de Avición CM 

* Tecnología de Avionos do Combate 
+ Tecnología do Hutecpteros 

+ Tocnología do Moteros a Roncción 
+ Elomantos do Diseno y Simulación 


Técnicas Enorgéticas 
* Técnicas do Combustón 


+ Enorgla Solar y de Hidtogoro, 
+ Lasor de Ala Enorgia 
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